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Abstract. Wc show that for a given bounded Apollonian circle packing 
V, there exists a constant c > such that the number of circles of 
curvature at most T is asymptotic to c- T a as T — > co. Here a ~ 1.30 is 

I ^ the residual dimension of the packing. For V integral, let 7r p (T) denote 

the number of circles with curvature bounded by T and prime. Similarly 
let (T) be the number of pairs of tangent circles with prime curvatures 
less than T. We obtain the upper bounds n v (T) < T a /logT and 

, , ^?(T) < T Q, /(logT) 2 , which are sharp up to constant multiple. The 

main ingredient of our proof is the equidistribution of expanding closed 

^\ horospheres in the unit tangent bundle of an infinite volume hyperbolic 

3-manifold r\H 3 under the assumption that the Hausdorff dimension of 
the limit set of V exceeds one. 
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Figure 1. Generations 1-3. 




Figure 2. A bounded Apollonian circle packing. 



1. Introduction 

A set of four mutually tangent circles in the plane with distinct points of 
tangency is called a Descartes configuration. Given a Descartes configura- 
tion, one can construct four new circles, each of which is tangent to three of 
the given ones. Continuing to repeatedly fill the interstices between mutu- 
ally tangent circles with further tangent circles, we arrive at an infinite circle 
packing. It is called an Apollonian circle packing, after the great geometer 
Apollonius of Perga (262-190 BC). 

See Figure 1 showing the first three generations of this procedure, where 
each circle is labeled with its curvature (that is, the reciprocal of its radius). 
Unlike the inner circles, the bounding circle is oriented so that its "outward" 
normal vector points into the packing. In Figure 2, the outermost circle has 
curvature —1 (the sign conveys its orientation). 

The astute reader would do well to peruse the lovely series of papers 
by Graham, Lagarias, Mallows, Wilks, and Yan on this beautiful topic, 
especially [16] and [15], as well as the recent letter of Sarnak to Lagarias 
[36] which inspired this paper. 
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Figure 3. An Apollonian circle packing bounded by two 
parallel lines. 

Counting circles in an Apollonian packing: Let V be either a bounded 
Apollonian circle packing or an unbounded one which is congruent to the 
packing in Figure 3. 

For V bounded, denote by N^(T) the number of circles in V in the 
packing whose curvature is at most T, i.e., whose radius is at least 1/T. For 
V congruent to the packing in Figure 3, one alters the definition of N^(T) 
to count circles in a fixed period. 

It is easy to see that N^{T) is finite for any given T > 0. The main 
goal of this paper is to obtain asymptotic formula for N V (T) as T tends to 
infinity. To describe our results, recall that the residual set of V is defined to 
be the subset of the plane remaining after the removal of all of the interiors 
of circles in V. Let a = a-p denote the Hausdorff dimension of the residual 
set of V . As any Apollonian packing can be moved to any other by a Mobius 
transformation, it follows that a does not depend on V . The current record 
towards the exact value of a = 1.30568(8) is due to McMullen [29]. 

Boyd [6] showed in 1982 that 

logA^(r) 
lim — — — = a. 

T^oo log T 

This confirmed Wilker's prediction [43] based on computer experiments. 

Regarding an asymptotic formula for N^(T), it was not clear from the 
literature whether one should conjecture a strictly polynomial growth rate. 
In fact, Boyd's numerical experiments led him to wonder whether "perhaps a 
relationship such as N V (T) ~ c-T a (log(T /c')) 13 might be more appropriate" 
(see [6] page 250). 

In this paper, we show that N V {T) has purely polynomial asymptotic 
growth. By f{T) ~ g(T) with T — > oo, we mean that lim-r^oc, = 1 

Theorem 1.1. Given an Apollonian circle packing V which is either bounded 
or congruent to Figure 3, there exists c = c(V) > such that as T — > oo, 

N V (T) ~ c-T a . 

In [6], Boyd actually considered the more general problem of counting 
those circles in a packing which are contained in a curvilinear triangle 1Z; let 
N^(T) count the circles having curvature at most T (see Fig. 4). For this 
question, it does not matter whether or not the full packing is bounded; the 
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Figure 4. Circles in a curvilinear triangle. 

counting function N n {T) is always well-defined. Since two such triangles 
are bi-Lipschitz equivalent, it follows from Theorem 1.1 that there exist 
constants ci , C2 > such that for all T S> 1, 

Cl -T a < N n (T) <c 2 -T a . 

Though we believe that the asymptotic formula N^(T) ~ c • T a always 
holds, our techniques cannot yet establish this in full generality. 

Primes and twin primes in an integral packing: A quadruple (a, b, c, d) 
of the curvatures of four circles in a Descartes configuration is called a 
Descartes quadruple. The Descartes circle theorem (see e.g. [9]) states that 
any Descartes quadruple (a, b, c, d) satisfies the quadratic equation 1 : 

(1.2) a 2 + b 2 + c 2 + d 2 = ^(a + b + c + d) 2 . 

Given any three mutually tangent circles with distinct points of tangency 
and curvatures a, b and c, there are exactly two circles which are tangent to 
all of the given ones, having curvatures d and d', say. It easily follows from 

(1.2) that 

(1.3) d + d' = 2{a + b + c). 

In particular, this shows that if a Descartes quadruple (a, b, c, d) correspond- 
ing to the initial four circles in the packing V is integral, then every circle in 
V also has integral curvature, as first observed by Soddy in 1937 [38]. Such 
a packing is called integral. 

It is natural to inquire about the Diophantine properties of an integral 
Apollonian packing, such as how many circles in V have prime curvatures. 



Arguably the most elegant formulation of this theorem is the following excerpt from 
the poem "The Kiss Precise" by Nobel Laureate Sir Fredrick Soddy [39]: 

Four circles to the kissing conic. / The smaller arc the bender. / 
The bend is just the inverse of / The distance from the center. / 

Though their intrigue left Euclid dumb / There's now no need for rule of thumb. / 
Since zero bend's a dead straight line / And concave bends have minus sign, / 

The sum of the squares of all four bends / Is half the square of their sum. 
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By rescaling, we may assume that V is primitive, that is, the greatest com- 
mon divisor of the curvatures is one. We call a circle prime if its curvature 
is a prime number. A pair of prime circles which are tangent to each other 
will be called twin prime circles. It is easy to see that a primitive integral 
packing is either bounded or the one pictured in Fig. 3. 

For V bounded, denote by ir^ (T) the number of prime circles in V of 
curvature at most T, and by tt^(T) the number of twin prime circles in V 
of curvatures at most T. For V congruent to the packing in Figure 3, one 
alters the definition of ir v (T) and tt^(T) to count prime circles in a fixed 
period. 

Sarnak showed in [36] that there are infinitely many prime and twin prime 
circles in any primitive integral packing V, and that 

* V{T) >> rHn^" 

(logT) 1 

Using the recent results of Bourgain, Gamburd and Sarnak in [5] on the 
uniform spectral gap property of Zariski dense subgroups of SL2(Z[i]), to- 
gether with the Selberg's upper bound sieve, we prove: 

Theorem 1.4. Given a primitive integral Apollonian circle packing V, 

(1) * v m « i&; 

(2) *?CT) « (iSf- 
Remarks: 

(1) The number of pairs of tangent circles in V of curvatures at most 
T turn out to be equal to 3N-p(T) up to an additive constant (see 
Lemma 2.5). Therefore, in light of Theorem 1.1, the upper bounds 
in Theorem 1.4 are only off by a constant multiple from the expected 
asymptotics. 

(2) One can also use our methods to give lower bounds for almost primes 
in a packing. A circle in V is called R-almost prime if its curvature 
is the product of at most R primes. Similarly, a pair of tangent 
circles is called R-almost twin prime if both circles are .R-almost 
prime. Employing Brun's combinatorial sieve, our methods show 
the existence of i?i,i?2 > (unspecified) such that the number of 
i?i-almost prime circles in V whose curvature is at most T is x ^f, 2 
and that the number of pairs of i?2-almost twin prime circles whose 
curvatures are at most T is x JJ^Prp ^ 



2 By f(T) x g(T), we mean g(T) « f(T) « g(T). 

•^The infinitude and in fact Zariski density of i?-almost prime and twin prime circles 
is obtained in [5]. The difference between our statements and theirs is that we can count 
with an archimedean norm. Neither method can specify Ri or R2, due to the lack of an 
explicit spectral gap, see [3]. 
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(3) A suitably modified version of Conjecture 1.4 in [5], a generalization 
of Schinzel's hypothesis, implies that for some c, C2 > 0, 



and 7r 2 (T) ~ c 2 



logT ~ " z v ' * (logT) 2 

(see the discussion in [5, Ex D]). The constants c and c 2 are detailed 
in [13]. 

Orbital counting of a Kleinian group in a cone: By Descartes' theorem 
(1.2), any Descartes quadruple (a, b, c,d) lies on the cone 
Q(x) = 0, where Q denotes the quadratic form 

1 



Q(a, b, c, d) 



a 2 + b 2 + c 2 + d 2 



-(a + b + c + dy 



In light of (1.3), one can "flip" the quadruple (a,b,c,d) into (a,b,c,d') via 
multiplication by S4, where 
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Let A denote the so-called Apollonian group generated by these reflections, 
that is A = (Si, S2, S3, S4), and let Oq be the orthogonal group preserving 
Q. One can check that A < Oq(Z) and that Q has signature (3, 1). There- 
fore A is a Kleinian group; moreover A turns out to be of infinite index in 
Oq(Z). 

For a fixed packing "P, there is a labeling by the Apollonian group A 
of all the (unordered) Descartes quadruples in V. Moreover the counting 
problems for N^(T) and N^iT) for V bounded can be reduced to counting 
elements in the orbit £p • A 1 C M 4 of maximum norm at most T, where 
is the unique root quadruple of V (see Def. 2.2 and Lemma 2.5). For V 
congruent to Figure 3, the same reduction holds with £p given by (0, 0, c, c). 

We prove the following more general counting theorem: Let 1 : PSL2(C) — ► 
SOp(IR) be a real linear representation, where F is a real quadratic form in 
4 variables with signature (3, 1). Let T < PSL2(C) be a geometrically finite 
Kleinian group. The limit set A(r) of T is the set of accumulation points 
of T-orbits in the ideal boundary 9oo(IHI 3 ) of the hyperbolic space H 3 . We 
assume that the Hausdorff dimension <5r of A(r) is strictly bigger than one. 



Theorem 1.5. Let vq G M 

with a discrete orbit vqT C 



be a non-zero vector lying in the cone F = 
Then for any norm \\ • \\ on M 4 , as T —> 00, 

#{v £ v T : \\v\\ < T} ~ c-T Sr 
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where c > is explicitly given in Theorem 7.1. 

There are two main difficulties preventing existing counting methods from 
tackling the above asymptotic formula. The first is that F is not required to 
be a lattice in PSL/2(C) (recall that the Apollonian group A has infinite index 
in Oq(Z)), so Patterson-Sullivan theory enters in the spectral decomposition 
of the hyperbolic manifold r\H 3 . The second difficulty noted by Sarnak in 
[36] stems from the fact that the stabilizer of vq in F may not have enough 
unipotent elements; in the application to Apollonian packings, the stabilizer 
is indeed either finite or a rank one abelian subgroup, whereas the stabilizer 
of vq in the ambient group G is a compact extension of a rank two abelian 
subgroup. 

In the similar situation of an infinite area hyperbolic surface, that is when 
T < PSL2(M), the counting problem in a cone with respect to a Euclidean 
norm was solved in Kontorovich's thesis [22], under the assumption that the 
stabilizer of vq in F is co-compact in the stabilizer of vq in PSL2(M). The 
methods developed here are quite different. Our approach to the counting 
problem is via the equidistribution of expanding closed horospheres on hy- 
perbolic 3-manifolds (see the next subsection for more detailed discussion). 
Our proof of the equidistribution works for hyperbolic surfaces as well, and 
in particular solves the counting problem in [22] for any norm and without 
assumptions on the stabilizer. In [21], we apply the methods of this paper to 
the problem of thin orbits of Pythagorean triples having few prime factors. 

Equidistribution of expanding horospheres in hyperbolic 3-manifolds: 

Let r be a geometrically finite torsion- free discrete subgroup of PSL^C). 
The main ingredient of our proof of Theorem 1.5 is the equidistribution of 
expanding closed horospheres in a hyperbolic 3-manifold r\H 3 . 

The group G = PSL2(C) is the group of orientation preserving isometries 
of the hyperbolic space H 3 = {(xi,X2,y) : y > 0}. The invariant measure 
on H 3 for the action of G and the Laplace operator are given respectively 
by 

dx\dx 2 dy , 9 / d 2 d 2 d 2 \ d 

and /\ = -y 2 [ + + )+y- 



y 3 \dy 2 dx\ dx\) dy 

Set 

1 x^ 




(1.6) 


N = {n x := | 




A = {a y := ( 




K = {g G G : 




fe ie 




M = < , 



fj : x € C}, 



= /} and 
: 6 G M}. 

By the Iwasawa decomposition G = NAK, any element g £ G can be 
written uniquely as g = n x a y k with n x £ N, a y 6 A, and k £ K. Via the 
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FIGURE 5 . Equi distribution of expanding horospheres 



map 



n x a y (0,0, 1) = (x,y), 



the hyperbolic space H 3 and its unit tang ent bundle T 1 (H 3 ) can be identified 
with the quotients G/K and G/M respectively. 

Denoting by [u] the image of u £ G under the quotient map G — > G/M, 
the horospheres correspond to iV-orbits [u]N = [uN] in G/M; note this 
is well-defined as M normalizes N. For a closed horosphere r\r[u]iV, the 
translates T\T[u]Na y represent closed horospheres r\r[na J/ ]A r which are 
expanding as y — > 0. See Figure 5. 

In the case of a hyperbolic surface r\H 2 of finite area, it is a theorem of 
Sarnak's [37] that such long horocycle flows are equidistributed with respect 
to the Haar measure. The equidistribution of expanding horospheres for 
any finite volume hyperbolic manifold can also be proved using the mixing 
of geodesic flows. This approach appears already in Margulis's 1970 thesis 
[25]; see also [12]. 

Alternatively, one can also use the measure classification for horospherical 
flows of Dani [11] and Ratner [33]. This has been carried out recently by 
Salehi Golsefidy and Mohammadi [30] for any finite volume homogeneous 
space of an arbitrary semisimple Lie group. 

In what follows we describe our equidistribution result for expanding 
closed horospheres on any geometrically finite hyperbolic 3 manifolds. 

Assuming 5r > 1, Sullivan [41] showed, generalizing the work of Patterson 
[31], that there exists a unique positive eigenfunction 0q of the- Laplacian 
operator A on r\H 3 of lowest eigenvalue 5r(2 — <5p) and of unit L 2 -norm, 
that is, J r \^ H 3 <Po{x, y) 2 ^dxdy = 1. Moreover the base eigen-value <5r(2 — 5v) 
is isolated in the L 2 -spectrum of A by Lax-Phillips [24] . 

Note that the closed orbit T\riV inside T^ryiHI 3 ) is an embedding of one 
of the following: a complex plane, a cylinder, or a torus. As 0o > 0, it is a 
priori not clear whether the integral 



Jn x £{Nnr)\N 

converges. We show that for any y > 0, the integral (a y ) does converge 
and is of the form 





00 K) = c 4>oV 2 &T + d <t>oV' 
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for some constants c^ > and d^ > 0. By f(y) ~ g(y) with y — > 0, we 

y^o g{y) 



mean that lim^o = 1- The measure dn on N is given by dn x = 



Theorem 1.7. Let T < G be a geometrically finite torsion-free discrete 
subgroup with 5r > 1 and let r\IW be closed. There exists e > such that 
for any ip £ C™(T\G) K = C C (T\H 3 ) ; and for all small y>0 



L 



ip(T\Tna y ) dn = C(j)0 ■ (ip, 0o)L2(r\H3) ■ y 2 5r (l + 0{y 6 )) 

(NnT)\N 

where the implied constant depends only on the Sobolev norm ofip. 

Thus as y — > 0, the integral of any function tp G C c (r\G)^" along the 
orthogonal translate T\TNa y converges to with the speed of order y 2 ~ 5r . 
It also follows that for ip £ C C (T\G) K , as y -> 0, 



tp(T\Tna y ) dn ~ (ip,<f>o) ■ / (poiTXTnay) dn. 

J (NnT)\N 



(Nr\T)\N J (NPiT)\N 

We denote by Or the set of vectors (p, v) in the unit tangent bundle 
T 1 (H 3 ) such that the end point of the geodesic ray tangent to v belongs to 
the limit set A(T) and by f2r its image under the projection of T 1 (H 3 ) to 

T 1 (r\H 3 ). 

Roblin [35], generalizing the work of Burger [7], showed that, up to a 
constant multiple, there exists a unique Radon measure fi on T 1 (r\H 3 ) 
invariant for the horosperical foliations which is supported on Cl-p and gives 
zero measure to all closed horospherical flows. 

In the appendix A written jointly by Shah and the second named au- 
thor, the following theorem is deduced from Theorem 1.7, based on the 
aforementioned measure classification of Burger and Roblin. In view of the 
isomorphism T 1 (r\H 3 ) = T\G/M, the following theorem shows that the 
orthogonal translations of closed horospheres in the expanding direction are 
equidistributed in T 1 (r\H 3 ) with respect to the Burger- Roblin measure fi. 

Theorem 1.8. For any ip £ C C (T\G) M , as y -> 0, 

/ ip(T\Tna y ) dn ~ c^ ■ y 2 ~ 5r ■ fi(ip) 

J(NnT)\N 

where fi is normalized so that £i((po) = 1- 

We conclude the introduction by giving a brief outline of the proof of 
Theorem 1.7. By Dal'bo (Theorem 3.3), we have the following classification 
of closed horospheres in terms of its base point in the boundary <9oo(IHI 3 ): 
r\riV is closed if and only if either oo ^ A(T) or oo is a bounded parabolic 
fixed point (see Def. 3.1). This classification is used repeatedly in our 
analysis establishing the following facts: 
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(1) For any bounded subset B C (N PI F)\N which properly covers 
(iV n r)\(A(r) - {oo}) (cf. Def. 4.5), 

<f>o(na y ) dn = 4> (na y ) dn + 0(y Sr ). 
'(Nnr)\N jb 

(2) Denoting by pb, € G C c (r\G) the e-approximation of B in the transver- 
sal direction, 



/ 

J Us 



[ M™*v) dn = Mo, PB,e) + 0{ey 2 ' 5 ^) + 0(/ r ). 

JB 



(3) For if; e C c (r\H 3 ), there exists a compact subset B = £>(supp(V0) C 
(iV n r) VV such that for all < y < 1 



(Arnr)vv 



if}(na y ) dn = if}(na y ) dn. 
Jb 



These facts allow us to focus on the integral of if; over a compact re- 
gion, say, B, of (N n r)\AT instead of the whole space. In approximating 
J B ip{na y ) dn with {a y ip, PB,e) , the usual argument based on the contracting 
property along the stable horospheres is not sufficient: the error terms com- 
bine with those coming from the spectral gap to overtake the main term! We 
develop a recursive argument which improves the error upon each iteration, 
and halts in finite time, once the main term is dominant. Using the spectral 
theory of L 2 (r\G) along with the assumption 5r > 1, we get a control of 
the main term of (a y ip,pB,e) as (ip,4>o) • {a y <j)o, pB,e)- 

For the application to Apollonian packings, we need to consider the max 
norm, which necessitates the extension of our argument to the unit tangent 
bundle, that is, the deduction of Theorem 1.8 from 1.7 as done in Appendix 
A. 

We finally remark that the power savings error term in (1) of Theorem 
1.7 is crucial to prove Theorem 1.4. 

Acknowledgments. We are grateful to Peter Sarnak for introducing us to 
this problem and for helpful discussions. We also thank Yves Benoist, Jeff 
Brock and Curt McMullen for useful conversations. 



2. Reduction to orbital counting 

2.1. Apollonian group. In a quadruple of mutually tangent circles the 
curvatures a, b, c, d satisfy the Descartes equation: 

a 2 + b 2 + c 2 + d 2 = -(a + b + c + d) 2 

as observed by Descartes in 1643 (see [9] for a proof). 

Any quadruple (a, b, c, d) satisfying this equation is called a Descartes 
quadruple. A set of four mutually tangent circles with disjoint interiors is 
called a Descartes configuration. 
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We denote by Q the Descartes quadratic form given by 

Q(a, b,c,d) = a 2 + b 2 + c 2 + d 2 -^(a + b + c + d) 2 . 

Hence v = (a, b, c, d) is a Descartes quadruple if and only if Q(v ) = 0. 
The orthogonal group corresponding to Q is given by 

Oq = {g G GL 4 : Q{v€?) = Q(v) for all v G M 4 }. 

One can easily check that the Apollonian group A := (Si, S2, S3, S4) defined 
in the introduction is a subgroup of Oq(Z) := Oq n GL^Z). 

Definition 2.1. (1) For "P bounded, denote by N^(T) the number of 
circles in V in the packing whose curvature is at most T, i.e., whose 
radius is at least 1/T. Denote by iV^T) the number of pairs of 
tangent circles in V of curvatures at most T. 

(2) For V congruent to the packing in Figure 3, N V (T) denotes the 
number of circles between two largest tangent circles including the 
lines and the largest circles. Similarly, N^iT) denotes the number of 
unordered pairs of tangent circles between two largest tangent circles 
including the pairs containing the lines and the largest circles. 

(3) For V bounded, denote by iP*{T) the number of prime circles in V 
of curvature at most T. Denote by ir^ (T) the number of twin prime 
circles in V of curvatures at most T. 

(4) For V congruent to the packing in Figure 3, one alters the definition 
of vr^'(T) and tt^(T) to count prime circles in a fixed period. 

We will interpret N^(T) and N£(T) as orbital counting functions on 
for a carefully chosen Descartes quadruple £ of V. 

Definition 2.2. A Descartes quadruple v = (a,b,c,d) with a + b + c + d > 
is a root quadruple ifa<0<6<c<(i and a + b + c > d. 

If V is bounded, Theorem 3.2 in [15] shows that V contains a unique 
Descartes root quadruple £ := (a, b, c, d) with a < 0. 

Theorem 2.3. [15, Thm 3.3] The set of curvatures occurring in V , counted 
with multiplicity, consists of the four entries in £, together with the largest 
entry in each vector £7* as 7 runs over all non-identity elements of the 
Apollonian group A. 

In [15], this is stated only for an integral Apollonian packing, but the same 
proof works for any bounded packing. Let be a non-returning walk 

away from the root quadruple £ along the Apollonian group, i.e., = 
^S^ ■ ■ ■ Sj n with Si k / Si k+1 , 1 < k < n — 1. Then the key observation in 
the proof of above theorem is that is obtained from li/"" 1 ) by changing 
one entry, and moreover the new entry inserted is always the largest entry 
in the new vector. 

This theorem yields that for T S> 1, 

N V (T) =#{ 7 Gi: llCf lima* < T} + 3. 
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Consider the repeated generations of V with initial 4 circles given by the 
root quadruple. Then a geometric version of Theorem 2.3 is that for n > 1, 
each reduced word 7 = Si n ■ ■ ■ Si 1 of length n corresponds to exactly one 
new circle, say C 7 , added at the n-th generation and the curvature of C 7 is 
the maximum among the entries of the quadruple £7* (cf. [16, Section 4]). 
Thus the correspondence (ft : 7 1— > C 7 establishes a bijection between the set 
of all non-identity elements of A and the set of all circles not in the zeroth 
generation, and the set {C 7 |7 7^ e, ||£7*|| < T} gives all circles (excluding 
those 4 initial circles) of curvature at most T. 

For each 7 7^ e in A, set 

<h{l) = {{C7, C 7 (l)}, {C 7 , C 7 (2)}, {C„ C 7 (3)}} 

where C 7 (i), i = 1, 2, 3, are the three circles corresponding to the quadruple 
£7* besides C 7 . Noting that each (C 7 ,C 7 (i)) gives a pair of tangent circles, 
we claim that every pair of tangent circles arise as one of the triples in the 
image of 4>2 , provided one of the circles in the pair does not come from the 
zeroth stage. If C and D form such a pair, they are from different stages of 
generation as no two circles in the same generation level touch each other, 
except for the initial stage. If, say, D is generated earlier than C, then for the 
element 7 £ A giving C = C 7 , which is necessarily a non-identity element, 
D must be one of C 7 (i)'s. This is because it is clear from the construction 
of the packing that every circle is tangent only to three circles from previous 
generations. 

Therefore </>2 yields a one to three correspondence from „4\ {e} to the set 
of all unordered pairs of tangent circles in V at least one circle of whose pair 
does not correspond to the root quadruple. Since there are 6 pairs arising 
from the initial 4 circles, we deduce that for T » 1, 

(2.4) JVf(T) = 3.#{ 7 Gi: Hfr'llma* < T} + 3. 

The above argument establishing (2.4) was kindly explained to us by Peter 
Sarnak. 

If V lies between two parallel lines, that is, congruent to Figure 3, there 
exists the unique c > such that V contains a Descartes quadruple £ := 
(0,0, c,c). 

In this case, the stabilizer of £ in A t is generated by two reflections S% 
and S\. One can directly verify that for all T S> 1, 

N V {T) = #{v G {A* : Hmax < T} + 3; 

and 

N%(T) = 3 • #{v € ZA* : Hmax < T} + 3. 

Lemma 2.5. Let V be either bounded or congruent to Fig. 3. 
(1) For allT^l, 

N V (T) = /# Stab ^(?) • #i v G £4* : Hmax < T} + 3 forV bounded 
I #{ v £ £A f : IMUax < T} + 3 otherwise. 
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(2) For allT^> 1, 

N%(T) = 3N V (T) - 6. 

(3) The orbit ^A l is discrete in M 4 . 

(4) For o«T» 1, 

4 

n V {T) <C = (wi,W2,W3,V4) S £„4* : ||u||max < F, is prime}. 

i=i 

(5) For o«T» 1, 

(T) <C ^ #{f G £4* : Ibllmax < T,Vi,Vj are primes}. 

Proof. The first two claims follow immediately from the discussion above, 
noting that the stabilizer of £ in .A* is finite for F bounded. The third claim 
follows from the fact that (T) < oo for any F > 0. For claims (4) and 
(5), note that for V bounded, 

ir v (T) <3 + #{ 7 £i: H^llmax is prime < T} 

4 

<C = (vi,v 2 ,v 3 ,v 4: ) £ £,A l : IMImax < T,Vi is prime}, 

i=i 

and 

^(T) < 5 + #{7 G .4 : ||^7*||max is prime < T, one more entry of £7* is prime} 
4 

< ^2^2#{ v = (v 1 ,V2,V3,v 4 ) £ £A l : IMImax < T,Vi,Vj are primes}. 

i=l j^i 

The claim (4) and (5) for V congruent to Fig. 3 can be shown similarly. 

□ 

We remark that there are bounded packings which are not multiples of 
integral packings: for instance, £ = (3 — 2y/3, 1, 1, 1) is a Descartes root 
quadruple which defines a bounded Apollonian packing. This is obvious 
from the viewpoint of geometry, but it is not at all clear a priori that the 
orbit ^A 1 should be discrete. Note also that there are other unbounded 
packings: by applying a suitably chosen Mobius transformation to a given 
packing, one can arrive at a packing which spreads uncontrollably to the 
entire plane, or one which is fenced off along one side by a single line. 

2.2. The residual set. We consider the upper half-space model for the 
hyperbolic space: 

M 3 = {( Xl ,x 2 ,y) eR 3 :y > 0} 
with metric given by dx i +d ^2+ d v _ j± <ii scre te subgroup of Isom(H 3 ) is called 
a Kleinian group. 

The ideal boundary d oc (M 3 ) of H 3 can be identified with the set of geodesic 
rays emanating from a fixed point xq 6 H 3 . The topology on ^(H 3 ) is 
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defined via the angles between corresponding rays: two geodesic rays are 
close if and only if the angle between the corresponding rays is small. 

In the upper half-space model, we can identify 9oo(BI 3 ) with the extended 
complex plane {(x±, X2,0)} U {00} = C U {00}, which is homeomorphic to 
the sphere § 2 . The space H 3 has the natural compactification H 3 U ^(H 3 ) 
(cf. [20, 3.2]). 

Definition 2.6. For a Kleinian group V, the limit set A(T) of F consists of 
limit points of an orbit Fz, z £ H 3 in the ideal boundary C = CU We 
denote by 5r the Hausdorff dimension of A(T). 

In this subsection, we realize the action of A on Descartes quadruples 
arising from V as the action of a subgroup, G^CP), of the Mobius trans- 
formations on the Riemann sphere C in a way that the residual set of V 
coincides with the limit set of G^(V). The residual set A(V) is defined to 
be the closure of all the circles in V, or equivalently, the complement in C 
of the interiors of all circles in the packing V (where the circles are oriented 
so that the interiors are disjoint). 

An oriented Descartes configuration is a Descartes configuration in which 
the orientations of the circles are compatible in the sense that either the 
interiors of all four oriented circles are disjoint or the interiors are disjoint 
when all the orientations are reversed. Given an ordered configuration T> 
of four oriented circles (Ci, C2, C3, C 4 ) with curvatures (61,62,63,64) with 
centers {(xj, yi) : 1 < i < 4}, set 



W v ■= 



(h 


61 


61X1 


612/1 \ 


62 


62 


62^2 


622/2 


63 


63 


63x3 


632/3 


\6 4 


64 


64X4 


642/4/ 



where 6j is the curvature of the circle which is the reflection of Ci through 
the unit circle centered at the origin, i.e., 6j = bi(xf + yf) — b~ l if 6, 7^ 0. 
If one of the circles, say Cj, is a line, we interpret the center (xi,yi) as the 
outward unit normal vector and set 6j = 6« = 0. 

Then by [16, Thm. 3.2], for any ordered and oriented Descartes configu- 
ration T>, the map g 1— > W^gW-p gives an isomorphism 

^:Oq(R)-Oq w (R) 
where Qw is the Wilker quadratic form: 



Qw = 



(0 


-4 





o\ 


-4 

















2 













V 



On the other hand, if M6b(2) denotes the group of Mobius transforma- 
tions and GM*(2) := M6b(2) x {±7} denotes the extended Mobius group, 
we have the following by Graham et al: 
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Figure 6. Action of A as Mobius transformations 

Theorem 2.7. [16, Thm. 7.2] There exists a unique isomorphism 

i:GM*(2)^0 0h ,(1) 
such that for any ordered and oriented Descartes configuration T>, 

(1) W lip) = Wv-k^)- 1 for any 7 G M6b(2). 

(2) W- V = -W v . 

Let Vq = (Ci, C2, C3, C4) denote the ordered and oriented Descartes con- 
figuration corresponding to the root quadruple £ of V. We obtain the fol- 
lowing isomorphism: 

$ Vo := 7T- 1 o : Oq(R) GM*(2). 

Denote by Sj := Sj(Po) the Mobius transformation given by the inversion 
in the circle, say, Ci, determined by the three intersection points of the circles 
Cj, j 7^ i. Figure 6 depicts the root quadruple (C\, . . . , C4) as solid-lined 
circles and the corresponding dual quadruple (Ci, . . . , C4) as dotted-lined 
circles. 

Note that Sj fixes Cj, j 7^ i and moves Ci to the unique other circle that 
is tangent to Cj's, j ^ i. 
We set 

G A (V) : = (Si, 32, 53,34). 
Lemma 2.8. For eac/i 1 < i < 4, 

$D (Si) =Sj; 

Zience 

$^(.4) = G A (V). 
Proof. If 7< := ^(Si), then by (2.7), 

On the other hand, by [16, 3.25], 

w s%[Vo) = s % w Vo . 

Therefore 

7l (P ) =s l (V ); 
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Figure 7. Fundamental domain 

and hence 7i = s$. □ 

Each inversion Sj extends uniquely to an isometry of the hyperbolic space 
H 3 , corresponding to the inversion with respect to the hemisphere whose 
boundary is C\. The intersection of the exteriors of these hemispheres is a 
fundamental domain for the action of G_a(V) on H 3 . 

The Hausdorff dimension of A(V), say a, is independent of V. Hirst 
showed in [18] that a is strictly between one and two. For our purpose, we 
only need to know that a > 1, though much more precise estimates were 
made by Boyd in [6] and McMullen [29]. 

Proposition 2.9. (1) Gjy(V) is geometrically finite and discrete. 
(2) We have A(G A (V)) = A(V), and hence a = S Ga{v) . 

Proof. For simplicity set T := G A {V). The group T is geometrically finite 
(that is, it admits a finite sided Dirichlet domain) by [20, Thm.13.1] and 
discrete since A is discrete and <£x> is a topological isomorphism. 

Clearly, T is non-elementary. It is well-known that A(T) is same as the 
set of all accumulation points in the orbit xq under V for any (fixed) Xq E C. 
On the other hand, by [16, Thm 4.2], A(7 3 ) is equal to the closure of all 
tangency points of circles in V and is invariant under T. This immediately 
yields 

A(r) cA(V). 

If xq G A(V), then any neighborhood, say U, of xq contains infinitely many 
circles. Since Ui<j<4r(Cj) is the set of all circles in V, there exist j and an 
infinite sequence 7$ G T such that ji(Cj) C U for all i. Therefore xo G A(r). 
This proves that A(V) C A(T). □ 

Set T-p to be the subgroup of holomorphic elements, that is, 

T v :=GU(P)nPSL 2 (C), 

since M6b(2) is the semidirect product of complex conjugation with the 
subgroup M6b + (2) = PSL^C) of orientation preserving transformations. 
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Then the above lemma holds with T-p in place of G^(V) as both properties 
are inherited by a subgroup of finite index. By the well known Selberg's 
lemma, we can further replace T-p a torsion-free subgroup of finite index. 

2.3. Reduction to orbital counting for a Kleinian group. Let G := 

PSL2(C) and r < G be a geometrically finite, torsion- free Kleinian sub- 
group. Suppose we are given a real linear representation i : G — ► SOf(IR) 
where F is a real quadratic form in 4 variables with signature (3, 1). As we 
prefer to work with a right action, we consider M 4 as the set of row vectors 
and the action is given by vg := i(g)v t for g G G and v G R . 

By Lemma 2.5 and the discussions in the previous subsection, Theorem 
1.1 follows from: 

Theorem 2.10. Suppose <5r > 1. Let vq G M 4 be a non-zero vector lying in 
the cone F = with a discrete orbit vqT C M 4 . Then for any norm \\ • \\ on 
M 4 , there exists c > such that 

#{v G v T : \\v\\ < T} ~c-T Sr . 

It also follows from Lemma 2.5 that this theorem implies 

N?(T) ~ (3c) • T a 

for the same c as in Theorem 1.1. 

3. Geometry of closed horospheres on T^ryH 3 ) 

Let G = PSL 2 (C) and H 3 = {(xi,X2,y) : y > 0}. We use the notations 
N, A, K, M, n x and a y as defined in (1.6). Throughout this section, we 
suppose that T < G is a torsion- free, non-elementary (i.e., A(T) consists of 
more than 2 points) and geometrically finite Kleinian group. 

Definition 3.1. (1) A point £ G A(T) is called a parabolic fixed point if 
£ is fixed by a parabolic isometry of T (that is, an element of trace 
±2). The rank of a parabolic fixed point £ is defined to be rank of 
the abelian group which stabilizes £. 

(2) A parabolic fixed point £ is called a bounded parabolic if it is of rank 2 
or if there exists a pair of two open disjoint discs in C whose union, 
say U, is precisely T^-invariant (that is, U is invariant by and 
7(CZ) n £/ / for all 7 G T not belonging to T ? ). 

(3) A point £ G A(r) is called a point of approximation if there exist a 
sequence 7» G T and z G H 3 such that 7«z —> £ and 7^21 is within a 
bounded distance from the geodesic ray ending at £. 

Since T is geometrically finite, it is known by the work of Beardon and 
Maskit [1] that any limit point £ G A(T) is either a point of approximation 
or a bounded parabolic fixed point. 

We denote by Vis the visual map from T 1 (H 3 ) to the ideal boundary 
d 00 (M 3 ) which maps the vector (p, v) to the end of the geodesic ray tangent 
to v. 



APOLLONIAN CIRCLE PACKING 



18 



Definition 3.2. • We denote by fV C T 1 (H 3 ) the set of vectors (p, v) 

whose the image under the visual map belongs to A(T). 
• We denote by fV the image of f2r under the projection T 1 (H 3 ) — > 

T 1 (r\H 3 ). 

We note that fir is a closed set, invariant under both geodesic flows 
and horospherical flows. For comparison, if V has finite co-volume, then 
A r = C U {00} and Q r = T 1 (r\H 3 ). 

For it 6 PSL2(C), we denote by [it] its image in PSL2(C)/M, which we 
may consider as a vector in T 1 (H 3 ). The horosphere determined by the 
vector [it] in r\T 1 (EI 3 ) corresponds to F\FuNM/M, which will be simply 
denoted by r\r[u]AT. We note that ix(oo) is identified with Vis [it]. 

Theorem 3.3 (Dal'Bo). [10] Let u £ PSL 2 (C). 

(1) If u(oo) £ A(T) is a point of approximation, then T\T[u]N is dense 
in fir- 

(2) The orbit T\TuN is closed in F\G if and only if either u(oo) ^ A(T) 
or it(oo) is a bounded parabolic fixed point for F. 

The first claim was proved in [10, Prop. C and Cor. 1] under the condition 
that the length spectrum of T is not discrete in M. As remarked there, 
this condition holds for non-elementary Kleinian groups by [17]. For the 
second claim, Dal'Bo only showed the equivalence with the closedness of the 
r\r[it]A r . However the above holds by the following lemma. 

Lemma 3.4. (1) We have F n NM = T f]N. 

(2) The orbit T\TNM is closed if and only if F\F N is closed. 

Proof If 7 G T n NM, but not in N, then Tr 2 (7) is real and < TV 2 (7) < 4 
where 1^(7) denotes the trace of 7. Hence 7 is an elliptic element. As F is 
discrete, 7 must be of finite order. However F is torsion-free, which forces 
7 = e. This proves the first claim. The second claim follows easily since the 
first claim implies the inclusion map Ff~)N\N — > rnA^M\A^M is proper. □ 

It follows that if r\rN is closed in r\PSL 2 (C), then r\rN is either 

(1) the embedding of a plane, if 00 ^ A(T); 

(2) the embedding of a cylinder, if 00 is a bounded parabolic fixed point 
of rank one; or 

(3) the embedding of a torus, if 00 is a bounded parabolic fixed point of 
rank two. 

For a subset X of the compactification H 3 u0 oo (BI 3 ), X denotes its closure 
in H 3 U5oo(M 3 ). 

Proposition 3.5. Assume that T\rA^ is closed. There exists a finite-sided 
fundamental polyhedron T C H 3 for the action ofF, and also a fundamental 
domain C C for the action of N n F such that for some r>l and for 
some finite subset Ir C F, 

{(xi,x 2 , y) e H 3 : x\ + ix 2 G Fn, x\ + x\ + y 2 > r} C U 7e / r 7JP". 
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Proof. Choose a finite-sided fundamental polyhedron T for T with oo G T. 
If oo ^ A(r), then oo lies in the interior of U^i^T for some finite Jr C 
r. As the exteriors of hemispheres form a basis of neighborhoods of oo in 
H 3 U <9oo(H 3 ), and T N = C, the claim follows. 

Now suppose that oo is a bounded parabolic fixed point of rank one. Let 
n Vl £ T fl # be a generator for m JV for v± G C, and fix a vector v 2 G C 
perpendicular to v±. The set J-pj := {s\Vi + s 2 w 2 G C : < s\ < 1} is 
a fundamental domain in C for the action of T n AT. By replacing T if 
necessary, we may assume that T C .Fat x M>o. Then by [26, Prop. A. 14 
in VI] , for all large c > 1 , the sets 

5(c) := {(sifi + s 2 V2,y) G H 3 : y > c or |s 2 | > c} 

is precisely invariant by TnA" and J-—S(c) is bounded away from 00. Since 
T n 5(c) / as 00 G T and 5(c) is precisely invariant by N n F, it follows 
that 

(3.6) {(xi,X2,y) G 5(c) : x\ + ix 2 G Jjv} C J 7 . 

Since 

{(xi,x 2 , y) : xi + ix 2 G F N , xj+xj+y 2 > r} C S(c) 

for some large c > 0, the claim follows in this case. 

If 00 is a bounded parabolic of rank two, S(c) := {(xi,x 2 ,y) : y > c} is 
precisely invariant by F n AT and F — 5(c) is bounded away from 00 by [26, 
Prop. A. 13 in VI]. Choose F N so that T C T N x M >0 . Then (3.6) holds 
for the same reason, and since Tn is bounded, 

{(x!,x 2 , y) : xi + ix 2 G J^v, + x 2 + y 2 > r} C 5(c) 

for some large c > 0. This proves the claim. □ 

Considering the action of A^ n F on 3oo(H 3 ) \ {00} = C, we denote by 
Ajv(r) the image of A(r) \ {00} in the quotient (N n T)\C, that is, 

A N (F) := {[n x ] G (N D r)\iV : x G A(r) \ {00}}. 

Proposition 3.7. IfF\FN is closed, the set Ajy(r) is bounded. 

Proof. This is clear if 00 is a bounded parabolic of rank two, as A^ n T\A^ is 
compact. If 00 ^ A(r), then A(r) is a compact subset of C, and hence the 
claim follows. In the case when 00 is a bounded parabolic of rank one, the 
claim follows from the well-known fact that A(r) lies in a strip of finite width 
(see [42, Pf. of Prop. 8.4.3]). This can also be deduced from Proposition 3.5 
using the fact that the intersection of the convex core of r\BI 3 and the thick 
part of the manifold r\HI 3 is compact for a geometrically finite group. □ 

Proposition 3.8. Assume that r\riV is closed. For any compact subset 
J C r\G, the following set is bounded: 

(3.9) A"(J) := {[n] G (N PI F)\N :F\FnAn J / 0}. 
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Proof. Let Tjsi and T be as in Proposition 3.5. If !Fn is bounded, there 
is nothing to prove. Hence by Theorem 3.3, we may assume that either 
oo ^ A(r) or oo is a bounded parabolic fixed point of rank one. 

To prove the proposition, suppose on the contrary that there exist se- 
quences rij £ F/v — ► oo, cij £ A, jj £ r and Wj £ J such that rijaj = ^jWj. 
As J is bounded, we may assume Wj — ► w £ G by passing to a subsequence. 
Then r yJ 1 rijaj — ► w/. Let 70 G T be such that 7ow(0,0,l) £ F. By the 
geometric finiteness of T, there exists a finite union, say J 7 ', of translates of 
F such that 7oWj(0, 0, 1) G F' for all large j. 

As rij — > 00, the Euclidean norm of njO,j(0, 0, 1) goes to infinity, and hence 
by Proposition 3.5, 

njdj(0, 0, 1) £ F' for all large j, 

by enlarging F' if necessary. 
Therefore for all large j, 

7077X^(0, 0, 1) = 70^(0, 0, 1) G F' n 7o77 1 (^0- 

Since J 7 ' Pi 7o7j~ 1 (F / ) 7^ for only finitely many 7/s, we conclude that {7?} 
must be a finite set. As rijaj = "YjWj £ 7jJ, njOj must be a bounded 
sequence, contradicting rij — > 00. □ 

Note that N(J) is defined so that for all y > 0, 
(iV n T)\Na y n J C iV(J). 

Corollary 3.10. Let ^ G C c (r\G) have support J. 

(1) TTie set N(J) defined in (3.9) is bounded. 

(2) For any function rj £ C C (N n r\iV) wii/i ?7|;v(j) = 1, we have for all 

y>o, 

/ ip(na y )dn = / ip(na y ) n(n) (in. 

J(Nnr)\N J(Nnr)\N 

Proof. The first claim is immediate from the above proposition. For (2), it 
suffices to note that ^(na y ) = for n outside of N(J), and hence 



ip(na y ) n(n) dn = ip(na y ) n(n) dn. 
Fy(./1 



I(N<1T)\N JN(J) 

Using n = 1 on N(J), the claim follows. □ 



4. The base eigenfunction o 

In this section, we assume that T < G = PSL2(C) is a geometrically finite 
torsion-free discrete subgroup and that the Hausdorff dimension 5r of the 
limit set A(r) is strictly bigger than one. Assume also that r\riV is closed. 

By Sullivan [41], there exists a positive F 2 -eigenfunction <fo, unique up to a 
scalar multiple, of the Laplace operator A on r\H 3 with smallest eigenvalue 
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<5r(2 — Sr) and which is square-integrable, that is, 

/ (f) (x,y) 2 \dxdy < oo. 
7r\e3 y 6 

In this section, we study the properties of 4>o along closed horospheres. 
For the base point o = (0, 0, 1), we denote by v a weak limit as s — ► 5r of 

— 3 

the family of measures on H : 

e L( , 70) E^ (0 ' T0) ^- 

The measure f G is indeed the unique weak limit of {v(s)} as s — > <5r- 
Sullivan [41] showed that f D is the unique finite measure (up to a constant 
multiple) supported on A(r) with the property: 

( 4-1 ) d(j*v ) = e _i ri 9 € ( 7 -i ,o) 

where /%(£i, £2) := liniz->£ 2) — d(z2, z) is the Busemann function and 
7*f G (A) := u (j(A)). The measure z/ , called the Patterson-Sullivan mea- 
sure, has no atoms. 

The base eigenfunction can be explicitly written as the integral of the 
Poisson kernel against u : 

(4.2) cj) G { x + iy)= f e- 5r ^ x+iy '°Uu (u) 

7ueA(r)\{oo} \ \\x-u\\ +V ) 
where A(r) is identified as a subset of C U {00}. The main goal of this 
section is to study the average of 4>o along the translates T\TNa y . 

Definition 4.4. For a given i/j 6 C(T\G) K , define 

ifj N (a y ) := / ip(na y ) dn. 

Jne(NnT)\N 

We will show that the integral 4>q (a y ) above converges, and moreover that 
all of the action in this integral takes place only over the following bounded 
set (see Proposition 3.7): 

Aiv(r) = {[n x ] g (N n r)\jv : x e A(r) \ {oo}}. 

If 00 is a bounded parabolic fixed point of rank one, then there exists 
a vector v\ G M? such that n vi is a generator for AfflT. Fixing vector 
i>2 perpendicular to v±, we arrive at a basis {^1,^2} for M 2 , i.e., we can 
decompose x G M? as x = xi^i + X2^2- 

Definition 4.5. We say that an open subset B C (NC\T)\N properly covers 
Ajv(r) if the following holds: 

(1) if 00 i A(r), then e (B) := iv£ ue A N (r),x<tB ||x - u\\ > 
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(2) if oo is a bounded parabolic of rank one, then 

€q(B) := inf \u 2 - x 2 \ > 0. 

xgB,ueA N (F) 

(3) if oo is a bounded parabolic of rank two, then B = (N n T)\N. 

Proposition 4.6. (1) (f>$(a y ) > y 2 ' Sr for all < y < 1. 

(2) For any open subset B C (iV n T)\A r which properly covers Ajy(r), 

(4.7) <${<*)= [ Mn x a y )dx + O eo{B) (y 5 n- 

Jb 

Proof. Choose a fundamental domain F/v C C for (NDT)\N. We first show 
that <f>Q (a y ) S> y 2_<5r for all < y <C 1. For a subset Q C Fat, the notation 
Q c denotes the complement of Q in Fat. Observe that for a subset Q C Fat, 



(4.8) / ^o(n x a y ) dx 
JxeQ 

= [ (||n|| 2 + l)^/ f- y—) T dxd v (u) 

= [ [I (\\uf + l) 5 r( - V dx) dv {u) 

J ueA(r) \JxeQ-u + lli/ll V / 

where the interchange of orders is justified since everything is nonnegative. 



X 



We observe that by changing the variable w 

(4.9) f ( V \ &T dx = y 2-Sv t ( 1 V r dw 

Jxm* VlMI 2 + ||y|| 2 / Jwm 2 \\\ w \\ 2 + 1 J 



y 2 ~ Sr 2ir [ — =- dr 

Jr>o (r 2 + l) ir 



<5 r - 1 

Case I: oo ^ A(r). In this case, we have 

loo = f (IMI 2 + l) 6v dv (u) < i/ (A(r)) < oo. 
j«GA(r) 

Therefore for Q = C, we obtain from (4.8) and (4.9) that 

(4.10) <j%{ay)= [ (H 2 + l)*&/ (u)- / (— — V r dx 
AeA(r) 7xeR 2 VlFH + 112/11/ 

= ^0 ' n ..2-Sr 

5 r -l 

This proves (1). To prove (2), let B be an open subset of C which properly 
covers Ajv(r). Then we have 

e := inf \\x - u\\ = mt\\\x\\ : x £ B c - A N (T)} > 0. 

ueA N (r),xeB*" 
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/ 4>o{n x a y ) dx <lj - ( ,, M9 - ,, M9 ) dx 

JxeBc J\\x\\>e \\\x\\ 2 + \\y\\ 2 J 

* " oy2 ~ 5r ■ I (ran) 



«5 r 

dw 



\\w\\>e y 
,2-6r 1 r 



27ruj y r • / r- dr 



aw 2 M '5 • i N 1 



y 2 



«/ r - 



Case II: oo is a bounded parabolic fixed point of rank one. 

We may assume without loss of generality that r Pi N is generated by 
(1, 0), that is, x i— > x + 1, and J 7 ^ = X2) '■ < x\ < 1,X2 G M}. There 
exists T > such that A(r) Clx [-T, T] by [42, Pf. of Prop. 8.4.3]. 

By computing the Busemann function, we deduce from (4.1) for k = 
(h,k 2 ) G rn N that 

/ lldl 2 + l \ Sr 
(4.11) d((n_ fc ).«/ )(u) = „ i - d«/ (u). 



_/c-u|| 2 + l^ 

We have by changing orders of integrations of xi and ui that 



(4.12) f Mx,y)dx 

Ju€Rx[-T,T]Jxe([0,l]xK.)-u\ \\ x \\ +V J 
= ! f ,, 2 ) r / {\\u\\ 2 + l)^du {u)dx. 

Jxm 2 \\\ x \\ +y J Jue\-x 1 ,l-x 1 }x\-T,T} 



Sr 

I 

E[-xi,l-xi] x[-T,T] 

We set ci and C2 to be, respectively, the infimum and the supremum of 

|2 , -i\<5r 



/ (||«f + 1)^^ («) 

Ae[-xi,l-xi]x[-T,T] 



over all x\ G R. We claim that 

< Ci < C2 < 00. 
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nui=ki+l— x\ 


rU2- 


=T 


u\=k\—x\ 


' U2 = 


-T 


ru\=k\-\-\—x\ 


I-U2- 


=T 


u\=k\—x\ 


1 112 = 


-T 



For x\ G [k\ — 1, k±), k\ G Z, by changing the variable u\ i— ► i*i + &i and 
recalling (4.11), we deduce 

/ (||«|| 2 + l)*-di/ («) 

Ju6[-xi,l-a;i]x[-T,T] 

/•ui=fci+l— />U2=T 

= / / (||n-fe|| 2 + l) 5 ^((n_ fc )^o)M 

J Ml =il - Si J uo = — T 

{\\u\\ 2 + l) 5 ^du {u). 

Since [fei — a;i, fci + 1 — xi] C [0, 2), we have 

ci > inf u ([k 1 -x 1 ,k 1 + l-x 1 ]x[-T,T]). 

fel6Z,xi6[fei — l,fci) 

Since the A?" n T-translates of [0, 1] x [— T, T] covers the support of v , we 
have u o ([0, 1] x [— T, T]) > 0. It then follows from the conformal property of 
v Q given by (4.11) that 

inf v ([ki — xi,ki + l — xi] x [— T,T]) > 

fei6Z,a;ie[fci— l,fei) 

and hence Ci > 0. On the other hand, 

(4.13) c 2 < sup i/ ([fcx -m,fei + 1 -xi] x [-T,T]) < oo 

fcieZ,a;i6[fci— l,fci) 

since f c is a finite measure. 

Therefore we deduce from (4.12) and (4.9) that 

/ (f) (x,y)dx^> ci / ( y — dx^>y 2 ' Sr , 
JxeT N JxeR 2 V IfII + y / 

proving (1). 

If i? is an open subset of ^jv = [0, 1] xl which properly covers Ajv(r), 
we have 

eo := inf I1X2 — X2I > 0. 
seB c ,«eA(r) 

Hence {(x, u) G M 2 x A^r) : 1 G 5 C - 11} is contained in the set 
{u\ G R, — 1/1 < xi < 1 — 1*1, \u2\ < T, |x2| > eo}. 

Since 



{1x1 G R, —iti < xi < 1 — 1*1} = {xi G R, — xi < 1*1 < 1 — xi}, 
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we deduce from (4.8) by changing the order of integrations that 
/ (j>o(x,y)dx 

_ f f (W±*) ir ^ o(u) 

AeA(r) JxeBc-u V \\ x \\ +V J 
< / H 2 / / (IM| 2 + l)* r cM^) dx 2 dxi 

■/xieR,|x 2 |>eo VIFII +2/ / \itii=-n Ju 2 =-T 



with C2 < oo as obtained in (4.13). 

The x\ integral can be evaluated explicitly and yields: 

<$ r -i/2 

J ^2 



A G B C r(<)r) V|xa|>eo V^2 +2/ / 

Li 



< v^-^r-l/2) , r /• /lV r - 1/2 fc 



« y 0T - 

Hence (2) is proved. 
Case III: oo is a bounded parabolic fixed point of rank two. In this case 
(2) holds for a trivial reason. But we still need to show (1). The argument 
is similar to the case II. Without loss of generality, we assume that N DT is 
generated by (1,0) and (0, 1), so that T N = [0, 1] x [0, 1]. 

Similarly to (4.12) we have 

<5 r 



<^K) = / f ll iia^ 2 ) " / (IMI 2 + V 5T dv {u)dx 

Jx&L 2 VIFII "+" V J Jtte[-Xl,l-Xl]x[-X2,l-X2] 

By (4.9), it suffices to show that 



c := inf / (||u|| 2 + l) &T du {u) > 0. 

xi,x 2 eM 1 / u6 [_ Xl) i_ Xl ]x[-x2,i-x2] 

For (xi,xz) S [ki — l,ki) x [k2 — 1,^ 2 ), by changing the variables Ui 
Ui + ki for ki G Z and by using (4.11), we deduce 

p«l = l — Xl /•U2 = l — 2:2 



rui=i— xi [-112=1-— X2 

/ / (||u|| 2 + l) fr <Mu) 



'ltl=— XI J U2=— X2 

/•Ul=fel + 1— Xl /•M2=fe2 + 1 — ^2 

= / / (|M| 2 + i)^ (u). 

Jui=ki—x\ JU2=k2—X2 

Observe that [&j — Xj, fcj + 1 — Xj] C [0, 2) and that f o ([0, l] 2 ) > 0, since the 
N n T-translates of [0, 1] x [0, 1] cover the support of v Q . Since ||n|| 2 + 1 is 
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bounded from below by a positive constant uniformly for all u G [0, l] 2 and 
ki G Z, we deduce that 

c > inf f ([A;i - Xj, fej + 1 - Xj]) > 0, 

again using the conformality of v Q . □ 

my y > 0, i/jere exist c<^ > 

0o Y K) = ^o2/ 2_5r +^o/ r - 



Corollary 4.14. For any y > 0, there exist c<f> > and a!^ > such that 



J/oo ^ A(r), i/ien d^ = 0. 

Proof. Since A0o = #r(2 — Sr)(j>o, it follows that 

As both and y 2 ~ Sr satisfy the above differential equation, we have 

^K) = c<t> y 2 ~ Sr + ^ / r 

for some c^ o ,d0 o > 0. Proposition 4.6 (1) implies that c<^ > 0. 

In the case when jF/v = C, the last claim is proved in (4.10). □ 

5. Spherical functions and spectral bounds 

We keep the notations set up in section 3. Let < s < 2, and consider 
the character Xs on the subgroup B := AMN of G defined by 

Xs{a y mn) = y s 

where a y G A is given as before, and m G M and n £ N. 

The unitarily induced representation (tt s := Ind^ Xs,Vs) admits a unique 
i^-invariant unit vector, say v s . 

By the theory of spherical functions, 



fs(g) ■= {n s (g)v s ,v s ) = v s {kg)dk 

Jk 



is the unique bi fT-invariant function of G with / s (e) = 1 and with Cf s = 
s(2 — s)f s where C is the Casimir operator of G. Moreover, there exist c s > 
and e > such that for all y small 

f s {a y ) = c s -y 2 - s {l + 0{tf)) 

by [14, 4.6]. 

Since the Casimir operator is equal to the Laplace operator A on K- 
invariant functions, this implies: 

Theorem 5.1. Let T < G be a discrete subgroup. Let 4> s G L 2 (T\G) K n 
C°°(T\G) satisfy A<fi s = s(2 — s)<f> 3 and \\4> s \\2 = 1- Then there exist c s > 
and e > such that for all small y < 1, 

(a y <f> s , s ) L 2 (r \ G) = c s ■ y 2 ~ s {l + 0{y 6 )). 
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In the unitary dual of G, the spherical part consists of the principal 
series and the complimentary series. We use the parametrization of s £ 
{1 + iR, [1,2]} so that s = 2 corresponds to the trivial representation and 
the vertical line 1 + iW corresponds to the tempered spectrum. Then the 
complimentary series is parametrized by V s , 1 < s < 2 defined before. 

Let {X{\ denote an orthonormal basis of the Lie algebra of K with respect 
to an ^4d-invariant scalar product, and define to := 1 — ^Xf. This is a 
differential operator in the center of the enveloping algebra of Lie(-fT) and 
acts as a scalar on each K- isotypic component of V s . 

Proposition 5.2. Fix 1 < sq < 2. Let (V,ir) be a representation of G 
which does not weakly contain any complementary series representation V s 
with parameter s > s$. Then for any e > 0, there exists c e > such that for 
any smooth vectors w\,w 2 £ V, and y < 1, 

\(a y w 1 ,w 2 )\ < c e -y 2 ~ s °~ e - |m«/i)|| • ||w(u> 2 )||- 
Proof. As a G-representation, tt has a Hilbert integral decomposition tt = 
IzeG®" 12 Pzdv(z) where G denotes the unitary dual of G, p z is irreducible 
and m z is the multiplicity of p z , and v is the spectral measure on G. By the 
assumption on ir, p z is either tempered or isomorphic to ir s for 1 < s < sq. 
As 1 < 3 — so < 2, there exists the complementary series (V3_ S0 , 7T3_ S0 ). 
We claim that the tensor product p z <S> tt3_ S0 is a tempered representation. 
Recall that an irreducible representation of G is tempered if and only if there 
exists a dense subset of vectors whose matrix coefficients are L 2+e -integrable 
for any e > 0. If p z is tempered, p z ^3- SQ is clearly tempered. If p z is 
isomorphic to tt s for some 1 < s < sq, and v z denotes the spherical vector 
of p z of norm one, then 

(a y (v z ® V3^ SQ ), (v z ®V3- S0 )) = {a y v z ,v z ) ■ (a y V3- So ,V3- So ) 

<0(y 2 - s -y 2 - (3 -^)=0(y) 

as s < so. 

Since the Haar measure on G satisfies d{k\a y k2) x y~ 3 dk\dydk2 for all 
< y < 1, it follows that the matrix coefficient (a y (v z l SiVs- S0 ), (v z <g)-y3_ S0 )) 
is L 2+e -integrable for any e > 0. This implies that the matrix coefficient 
of any linear combinations of g(v z ® t>3_ S(J ), g G G is also L 2+e -integrable 
for any e > 0. As they form a dense subset of vectors in the representation 
p z (g) 7T3_ So , it follows that p z (g) vr3_ So is tempered. 

Since tt <8> 7T3_ So = J z (B mz (p z ® , K3- SQ )dv(z), we deduce that it 7T3_ So is 
tempered. 

We now claim that for any e > 0, there is a constant c e > such that any 
iT-finite unit vectors u>i and W2, we have 

(5.3) \(a yWl ,w 2 )\ < c e -y 2 ~ s °- £ -nV(dim(K^))- 

i 

Noting that the -fT-span of Wi <S> v 3_ So has the same dimension as the K-span 
of Wi, the temperedness of tt ^3~ So implies that for any e > 0, there exists 
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a constant c e > such that 

{a y (wi (g> v 3 - So ), (w 2 <8> v 3 - So )) = (a y w 1 ,w 2 ) • (a y v 3 - SQ , v 3 - So ) 

< c e ■ y 1+e ■ [J V(&MKwi)). 

i 

As (a y vs- so ,vs- so ) = c- y~ 1+s ° (1 + 0(y^)) for some e > 0, the claim (5.3) 
follows. Passing from the above bounds of (6.6) for K- finite vectors to those 
for smooth vectors has been detailed in [28, Pf. of Thm 6]. In particular, in 
the case of G = SL2(C), the above degree of Sobolev norm suffices. □ 

Definition 5.4. For a geometrically finite discrete subgroup r of G with 
<5r > 1, we fix 1 < sp < Sp so that there is no eigenvalue of A between 
sp(2 — sp) and the base eigenvalue Ao = Sp(2 — Sp) in L 2 (r\G). 

By the theorem of Lax-Phillips [24], the Laplace spectrum on L 2 (T\G) K 
has only finitely many eigenvalues outside the tempered spectrum. Therefore 
1 < < <5r exists. The maximum difference between Sp and sp will be 
referred to as the spectral gap for T. 

Let {Zi, • • • , Zq} denote an orthonormal basis of the Lie algebra of G. 
Let T < G be a discrete subgroup of G. For / G C°°(T\G) n L 2 (T\G), we 
consider the following Sobolev norm S m (f): 

S m {f) = max{||Z n • • • Z lm {f)h ■ 1 < ij < 6}- 

Corollary 5.5. Let V be a geometrically finite discrete subgroup of G with 
Sp > 1. Then for any G L 2 (r\G) n C°°(T\G) K , i\) 2 G C™(T\G) and 
0<y<l, 

(ay^ifo) = ^i,^){a y M 2 ) + 0{y 2 -^S 2 ^ l )-S 2 ^2)). 

Here 4>o G L 2 (T\G) K is the unique eigenfunction of A with eigenvalue <5r(2 — 
5r) with unit I? -norm. 

Proof. We have 

L 2 (r\G) = w 5r ev 

where W$ r is isomorphic to Vs r as a G-representation and V does not con- 
tain any complementary series V s with parameter s > sp. Write ipi = 
(tpi, (l)o) 4>o + ipi- Since 4>o is the unique if-invariant vector in W§ T up to a 
constant multiple, we have ipj^ G V K . Hence by Proposition 5.2, for any 
e > and y < 1, 

(dyipufo) = (^i,0o> (0^0,^2) + (a y ipj-,ip 2 ) 

= (^i,<h){ay<h,ih) +0{y 2 ~ s rS 2 {i> 2 ) -S 2 ^ 2 )) 

since S 2 (tp^) <C <S(^). 

□ 
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6. EQUIDISTRIBUTION OF EXPANDING CLOSED HOROSPHERES WITH 
RESPECT TO THE BURGER-ROBLIN MEASURE 

Let T < PSL2(C) be a geometrically finite Kleinian group with St > 1. 
Assume that r\IW is closed. Let ji denote the Haar measure onG = NAK 
given by 

dfi(n x a y k) = y~ 3 dxdydk 
where dk is the probability Haar measure on K. We normalize <po so that 

1 

/r\H3 V 3 
By Corollary 4.14, we have 



!) (x,y) 2 -2dxdy = 1. 



/ 

J n 



(f>o(x, y) dx = c^y r + d^y r 

n x e(Nnr)\N 

where c^ > and d^ > 0. 

In this section, we aim to prove the following theorem. 

Theorem 6.1. For any V G C™(T\G) K , 

4>(n x a y ) dx = <Y>, fo) ■ c 0o ■ y 2 - 5r (l + 0{y^~ s ^)) 



I 



>n x e(Nnr)\N 

where the implied constant depends only on the Sobolev norm f or 
^ > 2(5 r -sr) ' ^ e v °l ume of N(supp(ijj)) and the volume of an open subset 
of (N n F)\N which properly covers Ajy(r). 

Most of this section is devoted to a proof of Theorem 6.1. 

Definition 6.2. For a given V G C°°(T\G) K and 77 G C c ((JVnr)\JV), define 
the function I^(V') on G by 

l v (^)(ay) : = / ip{n x a y )ri{n x ) dx. 

Jn x e(Nnr)\N 

We denote by iV~ the strictly lower triangular subgroup of G: 



The product map iV x i x A - x M — ► G is a diffeomorphism at a 
neighborhood of e. Let 1/ be a smooth measure on AN~M such that tin <8> 
v(an~m) = dfx. 

Fix a left-invariant Riemannian metric d on G and denote by U e the ball 
of radius e about e in C 

Definition 6.3. • We fix a non-negative function 77 G C^°((NnT)\N) 

with 77 = 1 on a bounded open subset of which properly covers 

Ajv(r). 
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• Fix eo > so that for the eo-neighborhood U eo of e, the multiplication 
map 

suppfa) x (U eo n AN~M) -► supp(r/)(?7 eo n AN~M) C r\G 

is a bijection onto its image. 

• For each e < eo, let r e be a non-negative smooth function in AN~M 
whose support is contained in 

w t := {u e n A)(u eo n N~)(u eo n M) 

and J w r e di^ = 1. 

• We define the following function p Vj€ on T\G which is outside 
supp(rj)U eo and for g = n x an~m G supp(»7)(i7 eo PI A/V _ M), 

Pr,, e (g) ■= V(n x ) ® r e (an~m). 
Recall from Proposition 4.6 that 

<(a s/ )=F ; (0 o )K) + O(y^). 
Proposition 6.4. T^e /iaue for all small < e <C eo and /or aZZ < y < 1, 

<tfi(ay) = (ay<h,Pv,e)v»(r\G) + O v (e ■ y 2 ^) + 0(/ r ) 
where the first implied constant depends on the Lipschitz constant for rj. 
Proof. Let h = a yi) n x m G W e . Then for n G AT" and y > 0, we have 

n/ia y = na yyo n~ x m. 

As the product map ixiVxif^Gisa diffeomorphism and hence 
a bi-Lipschitz map at a neighborhood of e, there exists £ > such that 
the e-neighborhood of e in G is contained in the product A^N^K^ of £e- 
neighborhoods for all small e > 0. 

Therefore we may write 

n yx = a yi n Xl h £ A eyeo N eyeo K 

and hence 

= n(ayy y 1 n Xl a yyoyi )a y y y 1 kim = n(n Xl yy oyi )ayy oyi kim. 
As 4>o is AT-invariant and dn is AT-invariant, 

/ M nha v) • n{n)dn = / 4>Q{n{n Xiy y oVl )ayy oyi ) ■ rj(n)dn 
J Nnr\N JNnr\N 

M na yy yi)(v(n) + 0{e))dn 

Nnr\N 

as r](n) — rj(nn') = O ri (e) for all n £ N and n'EJVfl U e . By Corollary 4.14, 
we deduce 



L 
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/ 4>o(nha y ) ■ rj(n)dn 

JNnr\N 

= / M na yy yi)v(n)dn + O v (e(f)$ (a yyoyi )) 
JNnr\N 

= c^iyyoyi) 2 -^ + ^ (yy yi) 5r + O^ey 2 '^) 

= c 0o y 2 -^(i + 0(1 - (y yi) 2 - 5r )) + O v (ey 2 ' s n + 0(y^) 

= C4>0 y 2 ~^{l + O{e)) + O{y s n 

as |y - 1| = 0(e) and \yi - 1| = 0(ye). 
As J r e du(h) = 1, we deduce 



(a y (/>o, Pr),e) = / r e (h) \ / 4> (nha y )ri(n) dn dv{h) 

Jw € \JNnr\N J 

= c 4>0 y 2 - s ^l + O 11 (e)) + O(y s n 
= ^(a y ) + O v (ey 2 - s n + 0(y 5 n. 



□ 



The following lemma is standard in Sobolev theory. 

Lemma 6.5. For ip G C^°(T\G) K , there exists $ G C™(T\G) K such that 

(1) for all small e > and h G U e , 

M9) ~ i>(9h) | < e ■ ${g) for all g G r\G. 

(2) S m (ip) <C S m +i(i/}) for each m G N, where the implied constant 
depends only on supp(V>)- 

Proposition 6.6. Let ^ G C°°(r\G)^. Then for any < y < 1 and any 

small e > 0, 

I MVOK) - (oj/V', Pr/,e)| < (e + y) • l v (tp)(a y ). 

Proof. If an-m G W e = (f7 e n A)(U eo n N~)(U £o n M), then 

(an~m)a y = a y a(a y -in~ a y )m G a^Wj 

since a y -i contracts iV~ by conjugation as < y < 1. 

As ^ is M-invariant, for any /t = an~m G W e , there exists an /i' G 
(f7 e n A)(U yeo n TV - ) such that 

\ip(na y ) — ip(nha y )\ = \4>{na y ) — ip(na y h')\ <C 4>(na y )(e + y). 
Hence 

^(ncty) — / ij){nha y )r e {h)dv{h)\ <C ip(na y )(e + y). 



Therefore 



Proof of Theorem 6.1: Recall that 77, eo, p v ,e = r]®r e are as in Def 6.3. For 
simplicity, we set p t = p v ,e- Noting that r t is essentially an e-approximation 
only in the yl-direction, we obtain that S 2 (p e ) = O^e" 5 / 2 ). 

We may further assume that 77 = 1 on iV(supp(V>)) by Corollary 3.10 so 
that 

/ ip(na y ) dn = l v (tp)(a y ). 

J(NnT)\N 

By Proposition 4.6, we also have 

<K) = i^ )(%) + o(/ r )- 

Set p = 5/2. Fix I G N so that 

^ (2-^ r )(p + i) ^ 
(<5 r - sr) 

Setting V'olfl') := t/HsO) we define for 1 < i < £, inductively 

Ma) ■= 4>i-i(g) 

where ipi-i is given by Lemma 6.5. 

Applying Proposition 6.6 to each tpi, we obtain for < i < I — 1 

l v (M( a y) = (a y ipi, p e ) + 0((e + y) l v (il> i+1 )(a y )) 

and 

Ir?(^)(%) = (a y ?Pe,Pe} + O v ((e + y)S 1 (tpe)) 
where the implied constant in the O n notation depends on J n dn. 
Note that by Corollary 5.5, we have for each 1 < i < £, 

{a y ipi,p e ) = {ipi,(f)o){a y (f) ,p e ) + 0(y 2 - Sr S 2 (pe)S 2 (il>i)) 

= O((a y 4> , p t ) ■ \\M\2) + 0(y 2 ^S 2 ( Pe )S 2 (M) 

= O((a y <f> , Pe ) ■ S t {$)) + 0(y 2 -^e'PS 2e ^)) 

Hence for any y < e, 

e-i 

l v W{a y ) = (a y ij, p t ) + 0{(a y ij k , p e )(e + y) k ) + O v (S e+1 (4>)(e + y) e ) 

k=i 

= (a y ^,p t ) + 0{(ay<h,P*)eSiW)) + 0{eS 2 ^)y 2 ~ s ^ e~ p ) + 0{S u ($)J) 

= ty, <M • {a y <fo, Pe) + O((a^o, Pe)e) + 0(y 2 - Sr e-P) + 0(e e ) 

= (^j, fa) • <f>N(ay) + 0(y 5 n + 0(ey 2 ~ s n + O^^) + 0(e e ) 

by Proposition 6.4, where the implied constants depend on the Sobolev norm 
S 2 ((ip) and J rjdn. 



APOLLONIAN CIRCLE PACKING 32 



(a y i/j,pe) L 2, r \ G )\ <C (e + y) • / ip{na y )i]{n)dn. 

J(Nnr)\N 



□ 
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Equating the two error terms 0(ey 2 ~ r ) and 0(y 2 ~ Sr e~ p ) gives the choice 
e = y ( 5 r-s r )/(p+i). By the condition on £, we then have e e <C y 2_<5r+ f ( 5r ~ Sr \ 
Hence we deduce: 

/ ^{na y ) dn = I v (^)(a y ) = <V>, <fo) ■ ^K) • (1 + 0(y^- Sr )). 

J{Nnv)\N 

Note that the implied constant depends on the Sobolev norm °f V' an d 

the L 1 -norm f rf dn, which in turn depends only on the volumes AT(supp(^>)) 
and any open subset of (TV n T)\N which properly covers Ajv(r). 

Burger-Roblin measure fi: In identifying ^(H 3 ) with K/M, we may 
define the following measure fi on T 1 (H 3 ) = G/M by 

dft(n x a y k) = y~ 1 ~ <5r dxdydv (k) 

for n x a y k G NAK/M, where v Q is the Patterson-Sullivan measure sup- 
ported on A(r) from section 4. The measure fi induces a Radon measure on 
T 1 (r\H 3 ) via the canonical projection and satisfies 

(6.7) A(V) = (V,<Ao) 

for AT-invariant functions ip G C c (T 1 (r\EI 3 )). 

Generalizing Burger's result [7], Roblin [35, Thm 6.4] proved: 

Theorem 6.8. The measure fi is, up to constant multiple, the unique Radon 
measure on T 1 (r\H 3 ) invariant for the horo spherical foliations whose sup- 
port is &r and which gives zero measure to closed horospheres. 

In Appendix A, the following theorem is deduced from Theorem 6.1. 

Theorem 6.9. For tp G C c (T 1 (r\H 3 )), 

ip(n x a y ) dx ~ C0 O • p,(ip) ■ y 2 ~ 5v as y ->■ 0. 



/ 

J n 



n x d(NC\T)\N 



7. Orbital counting for a Kleinian group 

Let 1 : G = PSL2(C) — > SOf(R) be a real linear representation where F 
is a real quadratic form in 4 variables of signature (3, 1). Let V < G be a 
geometrically finite torsion- free discrete subgroup with Sr > 1. Let vq G M 4 
be a non-zero (row) vector with F{vq) = such that the orbit vqT is discrete 
in R 4 . 

Since the orthogonal group Op(lR) acts transitively on the cone given by 
F = 0, there exists G Op(M) such that the stabilizer of 5^0 ^ s ec l ua l to 
N~M where N~ is the strictly lower triangular subgroup. In fact, g^VQ is 
unique up to homothety. Set 

To = 9o lv 9o- 

As for is discrete, it follows that ro\roA"M is closed. Hence by Lemma 3.4, 
the orbit T\TgoN is closed, equivalently Tq\ToN is closed. 
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Denote by (j) G L 2 (r \M 3 ) the unique positive eigenfunction of A with 
eigenvalue o~r(2 — Sr) and of unit L 2 -norm Jp ^ H 3 4>o{ x i y) 2 y~ 3 dxdy = 1. By 
Corollary 4.14, we have 

where c^, > and d^ > 0. Recall that the Patterson-Sullivan measure v Q 
on K which is normalized so that 

««,.)=/ ( j^ )* <M„). 

Theorem 7.1. For any norm || • || on M 4 , we /iaue, as T — > oo, 

#{v G Vo r : ||v|| < T} ~ o-p 1 • c 0() • f f Wvoigok-'g^W-^duoik)) ■ T Sr . 

\JkeK J 

V II ' II * s goK g^ 1 -invariant, then there exists e > such that 

#{v G v T : || v || <T} = «Jf 1 • <fo(e) • C( ^ • KIP^ • T 5r (l + 0{T^)) 

where e depends only on the spectral gap 5r — sr and the implied constant 
depends only on Ajv(r). 

By replacing T with To = g^Tgo, we may assume henceforth that go = e, 
and thus the stabilizer of vq in G is NM. By Lemma 3.4, The stabilizer of 
vo in r is simply T n N. 

Note that N~M fixes Vq, and that the highest weight (5 of the (irreducible) 
representation of l is given by (3(a y ) = y^ 1 . It follows that i(a y )vQ = y~ l v\, 
and hence v$a y = y _1 vo. 

Set 

B T := {v G voG : ||v|| < T}. 
Define the following function on T\G: 

f t{q) ■= XB T {v jg)- 
7 e(Afnr)\r 

Since vqT is discrete, Fx(g) is well-defined and 

F T (e) = #{v G v T : \\v\\ < T}. 
We use the notation: for ip G C c (r\G), 

ip N (a y ) := / tp(na y ) dn. 

J(NnT)\N 

Lemma 7.2. For any tp G C C (T\G) and T > 0, 

(F T ,V) = / / ^(a y )y- 3 dydk 

JkeM\K Jy>T- 1 \\v„k\\ 

where tjj k (g) = J meM tp(gmk)dm. 
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g e(Nnr)\G 



Proof. Observe: 

(F T ,ip) = XB T (v <yg)ip(g)dg 
r \ G 7 e(ATnr)\r 

XB T {v Q g)ip{g)dg 

= / / ip(n x ma y k)y~ 3 dx dmdy dk 

JkeM\K J y:\\v Q a y k\\<T Jn x me (NnT)\NM 

= [ [ (7 f ^.a^im^y-Hyik 

JkeM\K i J />T- 1 ||i-ofe|| \Jn x e(Nnr)\N JmeM J 

-f I (( Mn^y-'iyik 

JkeM\K J yyT^WvokW \Jn x e(Nnr)\N J 

= [ [ ^(a y )y- 3 dydk. 



>keM\K JyyT-^WvokW 

Define a function £„ : K — ► M. by 

£ V0 (k) = \\v k\\- 5 r. 
For ip £ C C (T\G), the convolution £„ * ip is then given by 

Zv *1>(9)= I ^{gk)\\v k\\-^dk. 



□ 



Ik&K 

Corollary 7.3. For any £ C c (r\G), we have, as T — ► oo, 

(F T ,4,)~8 r 1 -c <Po -T 5 r.jl(t V0 *ip). 
For ip £ C^°(T\G) K and \\ ■ \\ K-invariant, 

{F T , V) = W,<fo) ■ ^r 1 ■ • T 5 ? ■ \\v \\-^(l + Op-^-^)) 

where sr is as Def. 5.4 and the implied constant depends only on the Sobolev 
norm ofip, Ajy(r) and N(supp(ip)). 

Proof. Note that for any ip £ C C (T\G) and k £ K, the function ipk, defined 
in Lemma 7.2, is M-invariant. Applying Theorem 6.9 to ipk, we obtain that 
as y — ► 0, 

(7.4) / ip k (n x a y )dx ~ c^ • y 2 ~ Sr ■ fi{ip k )- 

Jn x &(NnT)\N 
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Hence by applying Lemma 7.2, inserting the definition of ipki and evalu- 
ating the y-integral, we get 



(F T ,ip) ~ C(t>0 • / / y-^K^k) dydk 

JM\K Jy>T- 1 \\v„k\\ 

= S~ 1 • c 0o ■ T Sr / / fi{ip{gmk)) ■ \\v k\\- 5r dmdk 

J M\K JM 



' M\K J M 

= ^ 1 -c^ -r <5r -/i(^ *^), 

proving the first claim. 

Now suppose that both -ip and the norm || • || are K- invariant. 

As tpk = ip, by Theorem 6.1, we can replace (7.4) by an asymptotic 
formula with power savings error term: 

/ iP{n x a y )dx = c^y 2 - 5 m>,<p )(l + O(y 2 7^~ s r))) 

Jn x £(Nnr)\N 

and the implied constant depends on the Sobolev norm of tp and Ajy(r). 
On the other hand, 

tvo**l>= INir^-V', 

and hence 

£(&o * VO = \\vo\\~ Sr ■ WAo)- 

Therefore 

(F T ^) = C(t)0 [ y- 1 - 5r (^,0 o )(l + 0(y 2 r^- s ^)) dy 

= Sf 1 ■ C(f>0 ■ T 5r • ||f |r 5r (l + 0(T-f( 5r - Sr ))) 

where the implied constant depends only on the Sobolev norm of ip and the 
set Ajv(r). □ 

Theorem 7.5. As T — »■ oo, 



F T (e) ~ Jf 1 • c^ • Qf_ |Nfc _1 |r* r <M&)) • T 5r . 

// || • || is K-invariant, then for some e > (depending only on the spectral 
gap 5 r - s r ), 

F T (e) = <5f 1 • Me) ■ c 0o ■ ||«o||-* r ■ T 5r (l + 0(T~ £ )). 
where the implied constant depends only on Ajv(T). 

Proof. For all small e > 0, we choose a symmetric e- neighborhood U t of e 
in G, which injects to T\G, such that for all T 3> 1 and all < e <C 1, 

B T U e C 5(i+ £ )t and 5 (1 _ e)T C n ueUe B T u. 

For e > 0, let <p e G C^°(G) denote a non-negative function supported on 
C/ e and with f G <j) £ dg = 1. We lift e to T\G by averaging over V: 

7er 
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Then 

(7.6) (F{i-*)T, $e) < F T (e) < (F {1+e)T , $ £ ). 
Note that 

(7.7) (F {1±e)T , $e) ~ ^f 1 • c*, • (T(l ± e)) 5r • A(^ * 
Considering the function R VQ : G — > M given by 

JkoeK 

for 5 = n x a y k, we have 

A(£«o* $ e)=/ / <S>e(gko)£, Vo {ko)dv (ko)d[i(g) 

JgGU e Jk eK 

= / / / $e(n x a y kko)£ vo (ko)dv (k )y~ 1 ~ 5r dxdydk 

J n x J <Xy J k J ko 

<$> e {n x a y k) (y 2 ~ 5r J iv l) {k 1 k)dv {ko)^ y~ 3 dxdydk 

<S>e{g)Rv (g) dg 
On the other hand, as / <J> e dg = 1 

/ $e(g)Rv (g) dg = R Vo (e) + 0(sup |i^ (x) - i^ (e)|) 

= [ |bo^ 1 |r <5r ^o(A;o) + 0(e). 



F T {e) ~ <5f 1 • • / ||«ofco ^-^di/oCfco) ■ ^ 



As e > is arbitrary, we deduce that 

-i(5p 

If || • || is if-invariant, we may take both U e and <p e to be i^-invariant. 
Hence by Corollary 7.3, we may replace (7.7) by 

(F (1±e)T , $ £ ) = <5 r 1 • C4>0 • \\v \r 5F • (T(l ± e)) 5r • (0o, * e > ■ (1 + 0{T^ ^-^))) 

where the implied constant depends on S2i(& e ) = ^ei^e) for some I > 1 
depending on the spectral gap of T, and the sets Ajy(r) and N(U e ) := {[n] G 

(ivnr)\iV:rnAn[/ e /0}. 

Since 

<0o, $ e ) = tfo(e) + 0(e), 

we have 

<F (1±e)T) $ £ ) = <5f 1 • c 0o • l^oir^ • T 5 r • </, ( e ) + 0(eT s n + 0(^7*-^-^)) 
for some q > depending on S2e(4> € )- Therefore by setting e 1+q = T~f^ Sr ~ Sr \ 
F T (e) = 5- 1 • C( ^ • \\v \\- 5 r • T 5r • Me) ■ (1 + 0(T~ e ')) 
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for e' = 7(i +q ) (^r — sr), where the implied constant depends only on Ajy(r). 

□ 

8. The Selberg sieve and circles of prime curvature 

8.1. Selberg's sieve. We first recall the Selberg upper bound sieve. Let 
A denote the finite sequence of real non-negative numbers A = {a n }, and 
let P be a finite product of distinct primes. We are interested in an upper 
bound for the quantity 

S(A,P):= £ °»- 

(n,P)=l 

To estimate S(A, P) we need to know how A is distributed along certain 
arithmetic progressions. For d \ P, define 

Ad := {a n £ A : n = 0(d)} 
and set |A d | := X)n=o(d) a «- We record 

Theorem 8.1. [19, Theorem 6.4] Suppose that there exists a finite set S 
of primes such that P has no prime factor in S. Suppose that there exist 
X > 1 and a function g on square-free integers with < g(p) < 1 for p\P 
and g is multiplicative outside S (i.e., g(d\d2) = g(d\)g(d2) if d\ and cfo are 
square- free integers with no factors in S) such that for all d \ P square- free, 

(8.2) \A d \=g{d)X + r d {A). 

Let h be the multiplicative function on square- free integers (outside S) given 
by h(p) = iz^p)- Then for any D > 1, we have that 

S(A,P)<xl M<0] + E r 3 (d)-\r d (A)\, 

\d<y/D,d\P ) d<D,d\P 

where Ts(d) denotes the number of representations of d as the product of 
three natural numbers. 

8.2. Executing the sieve. Recall from section 2 that Q denotes the Descartes 
quadratic form and A denotes the Apollonian group in Oq(Z). Fix a prim- 
itive integral Apollonian packing V with its root quadruple £ G Z 4 . Then V 

is either bounded or given by the one in Figure 3. 

To execute the sieve, it is important to work with a simply connected 
group. Hence we will set T_4 to be the preimage of SOq(M)° n A in the 
spin double cover Sping(IR) of SOq(M). Recall that a denotes the Hausdorff 
dimension of the residual set of the packing V . As shown in section 2.5, a is 
equal to 5r A , the Hausdorff dimension of the limit set of r_4. As SphiQ(IR) 
is isomorphic to SL^C), we have a real linear representation l : SL2(C) — > 
SOq(M) which factors through the quotient map SL^C) — > PSL^C). By 
setting T to be the preimage of T_4 under t, the counting results in the 
previous section are all valid for T. 
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As before, Bt denotes the ball in the cone: 

B T := {wGi 4 : Q{v l ) = 0, ||v|| < T} 

Since we are only looking for an upper bound, we may assume || • || is g^Kg^ 1 - 
invariant, where go 6 Oq(M) is such that the stabilizer of g^v^ is equal to 
N~M. We fix a small e > and let <j> e G C C °°(]H 3 ) = C™(G) K be as in the 
proof of Theorem 7.5 with G = SL 2 (C) and K = SU(2). 

Definition 8.3 (Weight). We define the smoothed weight to each 7 G V , 



w T {i)-= XB T {£ig)4>e{g)dn(g). 

Jg/k 

Let / be a primitive integral polynomial in 4 variables. Consider the 
sequence A(T) = {a n (T)} where 



a n {T) := Y, w ^)- 



7 estab r (5)\r 



Clearly, 



and 



|A(T)| = J> n (T)= Y W H7) 

n 7GStab r (C)\r 



(8.4) |A d (T)| = Y a n (T)= £ ^(7). 

For any subgroup To of T with 

Stab r (6 = Stabr (6, 

we define 

Fp(g) := Y XB T (^g), 
7 6Stabr(0\r 

and for each element 71 G T, we also define a function on ro\IH 3 by 

7er 

which is an e-approximation to the identity about [7f in ro\IH 3 . The 
function is simply the lift of </> e to ro\HI 3 and will be denoted by <I>f . 
For deZ, let r^(d) be the subgroup of T which stabilizes £ mod <i, i.e., 

r f (d) := {7 G T : £7 ^£(d)}. 

Note that 

Stab r (0 = Stabr e( d)(0- 
Lemma 8.5. (1) |A(T)| = (F$, ^> L2(r \H3) 5 
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(2) for any integer d, 

/«7l)=0(d) 

Proof. We have 

|A(T)| = ^ ™ T ( 7 ) 
7 estab r (0\r 



= E / XB T {iig)(i)e{g)dii{g) 



FHg)Ug) Mg) 

'G/K 

= [ F^(g)^(g) d»(g) 
Jy\g/k 

= { f ti $r }L2( r \ H 3) 

Expand (8.4) as 



IM^)I = E E «T(77i)- 



7ier £ (d)\r 7 eStab r (£)\r e (d) 

/«7l)=0(d) 



The inner sum is 

E «t(77i) E / Xs T (^7ff)0e(7r 1 5') 

76Stab r (€)\r«(d) 7 eStab r (0\r£(d) G/X 

JG/K 

Jrt:(d)\G/K 

Thus 

^ v ' \ 1 ,n /z,2 (r ((i )\H3) 

7 estabr(0\r e (d) v ev A 

□ 



Denote by Spec(r\H 3 ) the spectrum of the Laplace operator on L 2 (r\H 3 ). 
As mentioned before, the work of Sullivan, generalizing Patterson's, implies 
that Spec(r\H 3 )n[0, 1) / in which case Ao = a(2— a) is the base eigenvalue 
of A. For the principal congruence subgroup 

r(d) := {7 G r : 7 = I (mod d)} 

of T of level d, we have 

Spec(r\H 3 ) C Spec(r(d)\H 3 ). 
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The following is obtained by Bourgain, Gamburd and Sarnak in [5] and 
[4]- 

Theorem 8.6. Let T be a Zariski dense subgroup o/Sping(Z) with <5r > 1. 
Then there exist 1 < 9 < 5r such that we have for all square-free integer d, 

Spec(r(d)\M 3 ) n [6(2 - 9),5 r {2 - 5 r )} = {5 r {2 - S r )}. 

In order to control the error term for |A^(T)|, we need a version of Corol- 
lary 7.3 uniform over all congruence subgroups T^(d) of V. 

Proposition 8.7. There exists eo > 0, uniform over all square-free integer 
d, such that for any 71 6 T and for any congruence subgroup T^(d) ofT, we 
have 

' e ' 71 /L H r ddm s)- S r .[T:r^d)} m +U[€ ' 
for some d t > where the implied constant depends only on Ajv(T). 

Proof. Note that the congruence subgroup T(d) of level d is a finite index 
subgroup ofT^(d). Since 

Spec(r 5 (ti)\H 3 ) C Spec(r(ci)\H 3 ) 

the spectral gap Theorem 8.6 holds for the family T^(d), d square-free, as 
well. 

As we are assuming || • || is g^Rg^ 1 -invariant, by Corollary 7.3, we have 

(8.8) (f t ,*h*>) = (<^\^ {d v^; (d) -c^( d )-r^.|ieir-(i+o(T- £ °)), 

where eo depends only on the spectral gap for L 2 (T^(d)\E. 3 ) and the implied 
constant depends only on 

A N (g^(d)g Q ) = {[n x ] G (N D g l T^d)g )\N : x G A((fo ^(d)**))}- 
As r^(ci) is a subgroup of finite index in T, <5r 5 (cZ) = <^r and 

A(g 1 T^d)g )=A(g 1 Tg ), 
and moreover by the definition of T^(d), we have 

Stab r (0 = Stabr e( d)(0 > 

implying 

N nso^CtOso = NHg^Tgo. 

Hence 

An{9 ^{((fygo) = A N (gQ 1 Tg ), 
yielding that the implied constant in (8.8) is independent of d. 

By Theorem 8.6, eo can be taken to be uniform over all d. If To < T is a 
subgroup of finite index, the base eigenfunction in L 2 (ro\IH 3 ) with the unit 
L 2 -norm is given by 

<$> - 1 ^ 



v / [rTT ] 



'1— eo' 
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with cf)Q = Therefore 

This finishes the proof as 5r = a, by setting d e = (& e ,e, 4>o}- D 
Setting 

^ = j r : i . C0o . a ! e -iieir-T a , 

the following corollary is immediate from Lemma 8.5 and Proposition 8.7. 

Corollary 8.9. There exists eo > uniform over all square-free integer d 
such that 

where 

o° f (d)= £ i. 

7er e (d)\r 
/(£-y)=0(d) 

8.3. Proof of Theorem 1.4. We set 

fi(xi,X2,x 3 ,X4) = xi and /2O1, x 2 , x 3 , x 4 ) = £ix 2 . 
For <i square- free and i = 1,2, set 

5i (d) = o°(d)/[r : r e (d)]. 

Proposition 8.10. There exists a finite set S of primes such that: 

(1) for any square-free integer d = d\d2 with no prime factors in S and 
for each i = 1,2, 

9i(did 2 ) = gi{di) ■ gi(d 2 ); 

(2) for any prime p outside S, 

91 (p) G (0, 1) and 9l {p) = p- 1 + 0{p^' 2 ). 

(3) /or am/ prime p outside S, 

92(p) G (0, 1) and 53 (p) = 2^ + 0(jT 3 / 2 ). 

Proof. According to the theorem of Matthews, Vaserstein and Weisfeiler 
[27] , there exists a finite set of primes S so that 

• for all primes p outside S, T projects onto G(¥ p ) 

• for d = p\ ■ ■ ■ pt square-free with pi ^ S, the diagonal reduction 

r -► G{Z/dZ) -► G(¥ P1 ) x • • • G(¥ Pt ) 
is surjective. 
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Enlarge S so that G(¥ p )'s have no common composition factors for different 
p's outside S. This is possible because G = Spin(Q) can be realized as SL2 
over Q[\/— !]• Hence there exists a finite set S of primes such that for p 
outside S, 



(8.11) G(¥ p ) 



JSL 2 (F p ) xSL 2 (F p ) for p = 1(4) 
\SL 2 (F p2 ) for p = 3(4). 

It then follows from Goursat's lemma [23, p. 75] that T surjects onto 
G{Z/d{Z) x G{Z/d 2 Z) for any square- free d\ and d 2 with no prime fac- 
tors in S. This implies that for d = d±d 2 square- free and without any prime 
factors in S, the orbit of £ mod d, say O(d), is equal to 0{d\) x 0{d 2 ) 
in {Z/d{Z) 4 x (Z/d 2 Z) 4 = (Z/dZ) 4 . It also follows that O (d) is equal to 
0°(di) x O (d 2 ). Therefore g(d) = g(d 1 )g(d 2 ) as desired. 

Denote by V the cone defined by Q = minus the origin, i.e., 

4 4 

V = {(X!,x 2 , x 3 , x 4 ) + : 2 ^ ~ E x ^ = °^ 

i=i i=i 

Note that 

4 4 

Wi := {x G V : /i(x) = 0} = {(0, x 2 , x 3 , x 4 ) ± : ^ 2xf - (^ x,) 2 = 0}. 

i=2 i=2 

Since both quadratic forms 

44 44 
Q(xi,x 2 ,x 3 ,x 4 ) = ^2xf-(^Xj) 2 and Q(0, x 2 , x 3 , x 4 ) = ^ 2xf-(^ x { 

i=l i=l i=2 i=2 

are absolutely irreducible, we have by [2, Thm. 1.2.B], 

#V(F P ) = p 3 + 0(p 5 / 2 ), #Wi(F p ) = p 2 + 0(p 3 / 2 ). 

Since V is a homogeneous space of G with a connected stabilizer, by [32, 
Prop 3.22], 

0{p) = V(¥ p ), and hence O fl (p) = Wi(F„). 

Therefore we deduce gi(p) = #o\ p ) = p 1 + 0(p~ 3 / 2 ), proving (1). 

Now W 2 := {x G V : / 2 (x) = 0} is the union of two quadrics given by 
V n {xi = 0} and y n {x 2 = 0}. Hence 

#W 2 (¥ p ) = 2p 2 + 0(p 3 / 2 ). 

This yields that g(p) = 2p~ x + 0(p~ 3 / 2 ). □ 

We are now ready to prove Theorem 1.4. First consider /i(xi, x 2 , X3, x 4 ) = 
x x so that A(T) = {a n (T)} where 

a„(T) := Yl ^(7) 

7 sstab r (e)\r 

/l(?7)=« 
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is a smoothed count for the number of vectors (xi, X2, xs, x&) in the orbit 
of max norm bounded above by T and x\ = n. 
By Lemma 8.9 and as j^O^id) is multiplicative with #0 < j i (p) = p 2 + 

0(p 3 / 2 ), the quantity r^(A) in the decomposition (8.2) satisfies for some 
e > 0, 

r d (A) <^d 2 T a ~ e °. 

Thus for any E\ > 0, 

d<D,d|P 

which is < T°/logT for L> = T e °/ 4 , say. The key here is that D can be 
taken as large as a fixed power of T. Let P be the product of all primes 
p<D = T e °/ 4 outside of the bad set S. 

As h is a multiplicative function defined by h(p) = f° r P primes 

with <7i in Proposition 8.10, we deduce that (cf. [19, 6.6]) 

h(d) > logL> > logT, 

d<y/D,d\P 

and Theorem 8.1 gives 

S(A(T),P) «T Q /logT. 

Therefore 

#{(xi,x 2 ,x 3 ,x 4 ) e £A* : max < T, (x u \\ p) = 1} 

1<«<4 - LJ - 

p<T e /4 

<5(A((l + e)T),P) <T71ogT. 

Hence 

#{(^i, X2, X3, Xi) € ^A 1 : max < T, xi : prime} <C T a / logT. 

l<i<4 

Since this argument is symmetric in the Xj's, we have 

X2, xs, X4) G £.4* : max = prime at most T} <C T a /logT. 

l<i<4 

By Lemma 2.5, this proves 



* P < T » « log r- 



In order to prove 



(8-12) *?(T) « 

we proceed the same way with the polynomial /b^i, X2, X3, X4) = X1X2 and 
with the sequence A(T) = {a ra (T)} where 

a n (T) := Yl w T{l)- 

7 estab r («)\r 

/2(57)=™ 
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is a smoothed count for the number of vectors (xi, x 2 , xs, X4) in the orbit 
of max norm bounded above by T and x\x 2 = n. 
Note that g 2 (p) = 2p~ l + 0(p~ 3/2 ) by Proposition 8.10, and that 

/iW»(lo gj D) 2 »(logT) 2 

d<VD,d\P 

for h(p) = if.g^(p) f° r P primes. 
Therefore Theorem 8.1 gives 

S(A(T),P) <T a /(logT) 2 

which implies 

#{(xi,X2,x 3 ,X4) £ : max \xA <T,x±,x 2 : primes} <T a /(logT) 2 . 

l<j<4 

Again by the symmetric property of Xj's, we have 
#{x S £A* : max \xi\ < T, Xi, xj : primes for some j} < T a /(logT) 2 . 

1<«<4 

By Lemma 2.5, this proves 

^(T)«T7(logT) 2 . 
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A. Appendix: Non-accumulation of expanding closed 

HOROSPHERES ON SINGULAR TUBES (BY HEE Oh AND NlMISH SHAH) 

In this appendix, we deduce Theorem A.l from Theorem 6.1: Recall that 
r < G = PSL2(C) is a torsion free discrete geometrically finite group with 
Sr > 1 and that </>o G L 2 (r\BI 3 ) denotes the positive base eigenfunction of 
A of eigenvalue Sr(2 — 5r) and of norm f r \ m3 4>1)(g) dfj,(g) = 1. 

We continue the notations N,a y ,N~, etc., from section 6. We assume 
that r\riV is closed. By Corollary 4.14, for some > and > 0, 

/ (f) {na y )dn = c^y 2 '^ + d<p y Sr . 
Jr\TN 

We also recall the Burger-Roblin measure ft defined in Theorem 6.9 which 
is normalized so that /i(0o) = 1- 

Theorem A.l. For ip G C^T^ryH 3 )), 

ip(n x a y ) dx ~ 0^ ■ p,(ip) ■ y 2 ~ 5v as y ->■ 0. 



/ 

J n 



n x e(Nnr)\N 

Proof. The idea of proof is motivated by the proof of Lemma 2.1 in [34]. For 
each < y < 1, define the measure fi y on T\G/M = T 1 (r\H 3 ) by 

%M = c 0> _2+5r / ^{T\Y[n x a y }) dx 

J(NnT)\N 

for V G C C {T\G/M). 
Consider the family 

M := {Hy : < y < 1}. 

We claim that A4 is relatively compact in the set of locally finite Borel 
measures on T\G/M with respect to the weak*-topology. For any compact 
subset C C T\G/M, let if) be a if-invariant smooth non-negative function 
of compact support which is one over C. Then 

%(C) < nyty). 

As n y (ip) — ► A(V0 = (V 7 ) 00 ) by Theorem 6.1, the claim follows. 

It now suffices to show that every accumulation point of M. is equal to fi. 
Let /io be an accumulation point of M, which is clearly N- invariant. 

We denote by £p the set of vectors v G T 1 (r\H 3 ) the horospheres deter- 
mined by which is closed. In the identification of T 1 (r\H 3 ) with T\G/M, 
the set Sp corresponds to the image under the projection F\G — > T\G/M 
of the sets T\TgNA for T\TgN closed. 

Fix a bounded parabolic fixed point £o of T. If £o = <?o(°°) f° r 9o G 
PSL2(C), then a cusp, say, D(£o)j at £o is the image of U y>yi goNa y for 
some i/i > 1 under the projection ir : G — > T\G/M. 

There exists Co > 1 such that for any z = T\Tgona y G D(£o), 

(A.2) co 1 ^- 5 < Mz) < c y r - 5 
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where r G {1, 2} is the rank of £o (see [40, Sec. 5] as well as the proof of [8, 
Lem. 3.5]). 

Noting that Do = A)(£o) C £p, we first claim that for any weak limit, 
say, [i ,oiM, 

(A.3) MA)) = 0. 

Let Q be a relatively compact open subset of -Do- For y e > 1, setting 
-D e := Doa Vc C Do, we have that 



(A.4) f cgd^g) 

J D t K 



« y[- 2S . 



For the part of D e K inside the unit neighborhood of the convex core of V, 
this estimate follows from the proof of [8, Lem. 4.2]. When r = 1, the 
integral of <pQ over the part of D e K outside the unit neighborhood of the 
convex core of T is comparable to 



/*oo r 
J\ogy e Jx 



oo 

e~ 2St dxdt x y\- 2S . 

e* 



Fixing a neighborhood £7 of e in A" M such that DqU C -Do AT, the set 
D e Ua~^ is a neighborhood of Q. Therefore if \i Vi weakly converges to fio, 

c^ • Vo(4>o ■ XQ) < lim n I (<f>o ■ X DeUa -i)(na yi ) dn 



lim n ~2^5 I M na yi) ■ XD e u(na yi a ye ) dn. 



vt JV\TN 
As y e > 1 and U C N~M, 

Hence if na rae G D e U, then 

na^ G D e a y -i(a ye Ua y -i) C D e a y -iU C -D f7. 
Moreover if na We G TgQn'a y iMU for some n' G A and y' > yi, 
(A.5) na s G T g n a y > a y ^ (a y€ MU a y }) C g Q n'a yly -iMU. 

Using the formula 

0o(a) = / e-^^'°)^ (0 
^eA(r) 

for g G G (see (4.2)), it is easy to check that 

e -Sd(uo,o) Mg) < < e «5d( U o,o) 0o(5) 

for any g, u G G. Therefore we can deduce from (A. 2) and (A.5) that for 
some constant d > 1, 

4> (na y ) < c' ■ y 6 ~ r ■ 4>o(na y a ye ) 
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for all na y a Vc £ D e U. Hence 
C0 O • Mo (0o • Xq) 

< c 4>{) cy e 2 - r lim [ (^o ■ Xo e c/)(r\rna J/i a 2/e ) (in 

2/<^o (yiy e J z Jr\rA 

< c (j)0 c'y 2 ~ r / (0o ■ XD e K)(g) dn(g) by Theorem 6.1 

Jr\w3 

«y?" 2<5 by(A.4). 
Therefore 

MO(0O ■ XQ) < yl~ 2S - 

As <5 > 1, y e > 1 is arbitrary, and min g6 Q0o(<7) > 0, we have /J>o(Q) = 0. 
This proves the claim (A. 3). 

We now claim that Ho{£p) = for any weak limit hq: n yi — > mo- Suppose 
not. Since there are only finitely many cusps, there exist relatively compact 
open subset Q C £p and a bounded parabolic fixed point £o £ A(r) such 
that no{Q) > and its image a y (Q) = Qa y under the geodesic flow converges 
to £o as y — ► oo. Fix j/o > 1 such that 

c £> (£o)- 

By passing to a subsequence, we may assume n yiyo is convergent with a weak 
limit, say, fi' . Since Q^j/o c A)(£o)j by (A. 3), we have 

Mo(<2%o) = °- 

Therefore for any e > 0, there exists a neighborhood U e C N~M of e such 
that 

(A.6) // (Qa, ig<y 2+<5 e. 

Noting that Qa yo U t a~^ is a neighborhood of Q, we have 

0* < Km-lj j vxvN X QayoUt< {na yi )dn 

= VV (y^ / nrjv 

<e by (A.6). 
Since e > is arbitrary, Ho(Q) = 0. This proves 

Mo(£p) = 0. 
We deduce from Theorem 6.8 that 

for some qi > 0. On the other hand, by Theorem 6.1, 

Mo(0o) = A(0o) = 1- 
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It follows that a± = 1. □ 
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